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Abstract 

We describe our present understanding of the relations between the 
behaviour of asymptotically flat Cauchy data for Einstein's vacuum field 
equations near space-like infinity and the asymptotic behaviour of their 
evolution in time at null infinity. 
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1 Introduction 



There are no doubts any longer that the idea of gravitational radiation refers 
to a real physical phenomenon. Framing, however, a precise underlying math- 
ematical concept still poses problems. The work on gravitational radiation by 
Pirani, Trautman Sachs [S^j, [21] Bondi ^01) Newman and Penrose 
[35] and others, which was brought in a sense to a conclusion by Penrose |59| . 
|t)U|. is based on the idealization of an isolated self- gravitating system. It requires 
information on the long time evolution of gravitational fields which at the time 
could only be guessed. Ten years before these developments Y. Choquet-Bruhat 
had achieved a breakthrough in the mathematical analysis of the local Cauchy 
problem for Einstein's field equations However, the technical means to de- 
rive the fall-off behaviour of gravitational fields at far distances and late times 
from 'basic principles' were not available in the 1960's. In the meantime there 
has been a considerable progress in controlling the asymptotic structure of so- 
lutions to Einstein's field equations but it is still not quite clear which 'basic 
principles' to assume here. 

In the following we shall report on work which aims at closing various gaps in 
the study of gravitational radiation, the analysis of the Einstein equations, and 
the calculation of wave forms. Sections [5] -|5l present a fairly detailed discussion 
of the underlying analytical structures and of the recent results which led to the 
author's present understanding of the situation. To maintain the flow of the 
arguments, the reader is referred for derivations to the original literature. In 
sections El and will be given new results and detailed arguments. 

Penrose's proposal to characterize far fields of isolated systems in terms of 
their conformal structure ([Mj, [201) has been criticized over the years on sev- 
eral grounds; various variations, alternatives, etc. have been proposed (cf. [Ej, 
[T^. J7], [201, ESI J [10] I [Zn]i and references given therein). Some au- 

thors consider the smoothness requirements on the conformal boundary as too 
restrictive and suggest generalizations (cf. J7], [70], [75]). Doubts have been 
raised as to whether non-trivial asymptotically simple solutions to the vacuum 
field equations exist at all (Jl]) and it has been argued that the smoothness 
of the conformal boundary required in [55] excludes interesting physics ffO)'). 
The wide range of opinions on the subject is illustrated by the curious contrast 
between this emphasis on subtleties of the asymptotic smoothness and claims 
that 'null infinity is too far away for modelling real physics' (cf. j2S], |66| ) . 

In 120] even the asymptotically flat model is abandoned and replaced by a 
time-like cut model. The latter introduces a spatially compact time-like hyper- 
surface T which is chosen in an ad hoc fashion to cut off 'the system of interest' 
from the rest of the ambient universe. The idea then is to study the system 
which has thus been 'isolated' as an object of its own. 

The usefulness of any such suggestion can only be demonstrated by analysing 
its mathematical feasibility. This becomes clear when one tries to calculate 
wave forms numerically. Such calculations cannot be based on hand waving or 
physical intuition. The design of an effective numerical computer code requires 
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a precise mathematical formulation. 

The analysis of the time-like cut model reduces to a study of the initial 
boundary value problem for Einstein's field equations in which boundary data 
are prescribed on T and Cauchy data are given on a space-like hypersurface iS 
which intersects T in the space-like surface dS. In 021 has been given a fairly 
complete analysis of this problem for Einstein's vacuum field equations. This 
study is only local in time, but it provides insights into the basic problem. So 
far, the time- like cut model raises many more questions than it appears able to 
answer. 

How is T to be chosen ? Physical considerations may lead to suggestions 
when the system of interest is 'sufficiently far' away from other systems. How- 
ever, there is in general no preferred physical or geometrical choice for T. (It 
is instructive to compare this with the anti-de Sitter- type solutions, where the 
time-like boundary J" at space-like and null infinity is determined geometrically 
and the boundary data can be prescribed in covariant form (cf. ISH])-) 

The boundary must be characterized by some implicit or explicit geometrical 
condition. A natural choice is to prescribe its mean extrinsic curvature. Its 
evolution in time is then defined implicitly by a quasi-linear wave equation 
which itself depends in a non-local way on the data given on S and T (cf. |42|'). 
Long time calculations thus require an extra effort to control the regularity of 
the boundary. 

The gauge is related on the time-like boundary T directly to the evolution 
process. It depends on the (implicit) choice of a time-like unit vector field tan- 
gent to T. While the data which are prescribed on the space-like hypersurface 
S allow one to analyse the local geometry near S at any desired order, the data 
which can be prescribed on the boundary T provide very little information on 
the local geometry near T. All this makes it particularly difficult to show that 
the gauge and the constraints are preserved under the evolution in time. 

These properties imply in general a non-covariance of the boundary condi- 
tions and data. Moreover, due to the fact that no causal direction is distin- 
guished on T there does not seem to exist a natural 'no incoming radiation 
condition' and, in particular, no natural concept of 'outgoing radiation'. In 
fact, it appears difficult to associate with the initial boundary value problem 
any 'simple' quantities which characterize the system and its dynamics and 
which can be related to observational data. 

While the discussion in 021 singles out data which arc mathematically ad- 
missible, it is far from clear what should be prescribed on T from the physical 
point of view. The 'correct' data induced by the ambient universe will never be 
known. The information fed into 'the system' by the data prescribed on T can 
hardly be assessed. In long time calculations it may alter the character of the 
system drastically. 

Because of these difficulties the time-like cut model appears not very promis- 
ing. Nevertheless, it is of interest because of its similarity to the standard 
approach to numerical relativity, where an artificial time-like boundary is in- 
troduced to render the computational grid finite. It is expected here that the 
assumption of asymptotic flatness together with a judicious choice of the bound- 
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ary will alleviate some of the difficulties pointed out above. 

At present the only satisfactory solution to the gravitational radiation prob- 
lem is based on the assumption of asymptotical flatness and the most elegant 
and geometrically natural definition of the latter is provided by the idea of 
the conformal boundary at null infinity introduced in |59) . While useful physi- 
cal concepts can be associated with a conformal boundary which is sufficiently 
smooth (cf. 1^, US and the references given there), the possible degree of 
differentiability, which encodes the fall-off behaviour of the gravitational field, 
still poses questions. This article deals with this particular issue and tries to 
disentangle its various aspects and difficulties. 

Einstein's field equations admit certain conformal representations which in 
the following will be referred to as conformal field equations. These equations 
are 'regular' in the sense that they imply in a suitable gauge equations which 
are hyperbolic even at points of null infinity (|2H1, EH)- This fact has been 
used to show that the smoothness of the conformal boundary is preserved if it is 
guaranteed on the initial slice S of an hyperboloidal initial value problem f|31|. 
cf. also ISSI). The subsequent analysis of hyperboloidal initial data ([3], 
showed the existence of a large class of smooth hyperboloidal data for 
the conformal field equations. The construction of such data requires the 'free 
data' to satisfy a finite number of conditions at the space-like boundary dS at 
which the hyperboloidal slice S intersects future null infinity J''^. 

However, the work referred to above also shows the existence of a large class 
of hyperboloidal data which are smooth on 5 \ dS but possess a non-trivial 
polyhomogeneous expansion at 95, i.e. an asymptotic expansion in terms of 
rf.k jQgj ^ where a; is a defining function of the boundary dS^ which vanishes on 
dS. Logarithmic terms can occur as a consequence of the constraint equations 
even if the free data extend smoothly to dS. Recently, it has been shown 
that certain hyperboloidal data which are polyhomogeneous at dS evolve into 
solutions to the conformal field equations which possess generalized conformal 
boundaries near the initial slice ([TH|, |S7]). While the precise behaviour of 
these solutions near that boundary still needs to be analysed, the result shows 
that the use of the conformal field equations and the characterization of the 
edge of space-time in terms of its conformal structure are not restricted to 
asymptotically regular situations. 

We conclude from these results that in the standard Cauchy problem the 
field equations decide on the degree of smoothness of the conformal boundary 
at null infinity in arbitrarily small neighbourhoods of space-like infinity. 

There are other reasons to study the region near space-like infinity. The 
hyperboloidal initial value problem is intrinsically time-asymmetric. To anal- 
yse in the same picture incoming radiation, a non-linear scattering process, and 
outgoing raditation, one needs to include space-like infinity (as pointed out al- 
ready in (HOI)- Also, if the hyperboloidal data are not distinguished by special 
features as, for instance, the presence of a trapped surface, it is not clear which 
part of the imagined space-time is covered by their evolution. They could repre- 
sent a hypersurface close to time-like infinity or close to a Cauchy hypersurface 
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(a difficulty shared with the characteristic initial value problem and the initial 
boundary value problem). 

This should not obscure the fact that numerical calculations of space-times 
from hyperboloidal data allow one to determine wave forms for many 'realis- 
tic' physical processes. So far the only semi-global calculations of space-times, 
including their radiation fields at null infinity, are based on hyperboloidal and 
characteristic initial value problems (cf. |2ZI, [H], and the article by L. 
Lehner and O. Reula, this volume). 

We are thus left with the following task: (i) characterize the data which 
evolve near space-like infinity into solutions of prescribed smoothness at null 
infinity, (ii) analyse for which of these data physical concepts and requirements 
(linear and angular momentum at space-like and null infinity, reduction of the 
asymptotic gauge (BMS) group to a Poincare group, . . . ) can be meaningfully 
introduced and a satisfactory physical picture can be established. 

The first step is technically the most difficult one. It requires us to control 
under fairly general assumptions the effect of the quasi-linear, gauge hyperbolic 
field equations over infinite regions of space-time. Moreover, the asymptotic 
behaviour of the solutions has to be determined with a precision which excludes 
any further refinement. 

Once this step has been taken, many considerations of the second step will 
reduce to straightforward, though possibly quite lengthy, calculations. Nev- 
ertheless, the second step is of crucial importance. At this stage one has to 
observe that the notion of asymptotic flatness is not part of the general theory; 
it is an idealization which chosen to serve a purpose. While it is suggested to 
us by important solutions such as those of the Kerr family, it is far from being 
determined by the equations alone. There remains a large freedom to decide on 
the asymptotic behaviour of the fields. 

To make one's choice, one needs to know the mathematical options and 
has to decide on the physical questions to be answered. A theorem which 
characterizes the most general Cauchy data on 5 = R'^ for which the maximal 
globally hyperbolic Einstein development is null geodesically complete and for 
which the Riemann tensor goes to zero at (null) infinity would be mathematically 
quite an achievement but, by itself, insufficient from the point of view of physics. 

We are not interested here in discussions of asymptotically flat solutions 
with 'observations' referring to the roughness of the asymptotic structure as, for 
instance, in IJO]- We rather wish to understand whether (i) the solution models 
a 'system of physical interest' and (ii) its far field and asymptotic structure 
allow one to extract information on the system which characterizes its physical 
nature and can be related to observational data. 

This task is neither easy nor well defined. The studies of the last 40 years 
provide some understanding of the situations one may expect to observe (col- 
lapse to a black hole, mergers of black holes, . . .). By exploring, however, the 
questions above in a general setting, new phenomena may be encountered (cf. 
|12| . for an example). But given that the interior is understood to some 
degree, what do we do about (ii) ? 

Recent results on the constraint equations exhibit possibilities to modify 
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asymptotically flat vacuum data 'far out' without affecting the interior. The 
data can be made to agree near space-like infinity with exact Schwarzschild or 
Kerr data {^U\, |23)i with even more general static resp. stationary data, or 
with data which are only asymptotically static resp. stationary (^|) (cf. also 
the discussion in the article by R. Bartnik and J. Isenberg, this volume, for 
other techniques of modifying or extending solutions to the constraints) . 

These results have been used to settle a question which has been open for 
a long time. Since data which are static or stationary near space-like infinity 
evolve into solutions which possses a smooth conformal boundary at null infinity 
(cf. 123) J these solutions contain smooth hyperboloidal hypersurfaces. Recently 
P. Chrusciel and E. Delay have shown the existence of families of Cauchy data on 
R'^ which are static outside a fixed radius and have members of arbitrarily small 
ADM-mass. The corresponding solutions contain hyperboloidal hypersurfaces 
to which the results of apply- This demonstrates the existence of non- 
trivial asymptotically simple solutions to Einstein's vacuum field equations with 
prescribed smoothness of the asymptotic structure ([ISp. 

More recently S. Klainerman and F. Nicolo revisited their work in |5^) and 
showed (HSI) that their solutions will have the Sachs peeling property ((H3|, E^) 
if the data are subject to certain asymptotic conditions. However, the class of 
data which meet these requirements still needs further analysis. 

The new flexibility in constructing asymptotically flat initial data also allows 
one to illustrate some difficulties of the asymptotically flat space-time model. 
Let (5, d) denote the initial data where S is the hypersurface considered in our 
discussion of the time-like cut model and d indicates the fields induced on S by 
the cosmological model. Suppose that {S',d') is an asymptotically flat initial 
data set for which there exists an embedding (j) : S S' such that the push 
forward of d hy cj) is in a suitable sense 'close' to d' on (l){S). The evolution in 
time of the data (5', d') can then be considered in some neighbourhood of 4>{S) 
as a good approximation of the evolution of d in the cosmological model. 

If the set S is chosen large enough and close to the region where the system 
is undergoing a wave generation process, the main part of the wave signal will 
reach null infinity at a finite retarded time. The fact that for very late times the 
data on S' \ <j){S) will create a deviation of our solution from the cosmological 
one is likely to be irrelevant in many interesting situations. From a pragmatical 
point of view it may be considered the main purpose of the asymptotically 
flat space-time extension beyond the domain of dependence of (f>{S) to allow 
perturbations of the gravitational field generated near (j){S) to unfold into a 
clean wave signal which can be read off at null infinity. 

Since changes near space-like infinity affect the field, however weakly, at all 
later times, they may have an important effect in the case of black hole solutions. 
One may envisage the collapse of pure gravitational radiation to a black hole 
as being modelled by vacuum solutions which arise from smooth asymptotically 
flat data on M^, admit smooth, complete (cf. @7]) conformal boundaries J'^, 
and possess future event horizons while all past directed null geodesies require 
endpoints on J^^ . At present nothing is known about such solutions and they 
may not exist. Is it conceivable then that there exist solutions which show 
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all the (suitably generalized) features listed above but have a rough conformal 
boundary ? Could such boundaries allow for a 'higher radiation content' ? 
If that were the case the restriction to smooth conformal boundaries might 
preclude the discussion of certain interesting physical phenomena. 

Clearly, the large freedom in constructing asymptotically flat extensions 
should neither be used to import irrelevant information into the system nor 
to suppress important features. The replacement of an extension by one which 
is strictly Kerr (say) near space-like infinity introduces a transition zone on 
the initial hypersurface which mediates between the given and the Kerr data. 
The resulting 'wrinkles' in the solution are recorded in the radiation field at 
null infinity. Is this information physically insignificant or does it indicate that 
something important has been ironed out by forming the new extension ? 

This question points again to the need of understanding the detailed be- 
haviour of the fields near space-like infinity. In the standard conformal represen- 
tation space-like infinity with respect to the solution space-time is represented 
by a point, usually denoted by With respect to an asymptotically flat Cauchy 
hypersurface S space-like infinity is then also represented by a point, denoted 
by i, which becomes under conformal compactification a point in the extended 
3-manifold S = S U {i}. Unfortunately, if mAOM 7^ the conformal initial data 
are strongly singular at i. This is the basic stumbling block for analysing the 
field near space-like infinity in terms of the standard conformal rescaling. 

The constraint equations on the Cauchy hypersurface S impose only weak 
restrictions on the asymptotic behaviour of the data near space-like infinity. It 
is easy to construct data which at higher orders will become quite 'rough' near 
i and which can be expected to affect the smoothness of the fields near null 
infinity in a physically meaningless way. Thus one will have to make a reason- 
able choice and find a class of data which allows one to perform a sufficiently 
detailed analysis of their evolution in time while still being sufficiently general 
to recognize basic features of the asymptotic behaviour at null infinity. 

In the following it is assumed that the data on S represent a space-like slice 
of time reflection symmetry, so that the second fundamental form vanishes, and 
that their conformal structure, represented by a conformal 3-metric h on S, 
extends smoothly to the point i. Only these conditions will be used in the 
following discussion, no a priori assumptions on the evolution in time will be 
made. We note that the assumptions are made to simplify the calculations, 
there exists a large space for generalizations. 

Somewhat unexpectedly, the attempt to analyse for data as described above 
the evolution near space-like infinity i in the context of the conformal field 
equations led to a finite regularization of the singularity at space-like infinity 

m)- 

In a certain conformal scaling of the conformal initial data the choice of 
the frame and the coordinates is combined with a blow-up of the point i to 
a sphere X° such that the initial data and the gauge of the evolution system 
become smoothly extendable to in a different conformal scaling. Moreover, 
the general conformal field equations imply in that scaling a system of hyperbolic 
reduced equations which also extends smoothly to fscction [5.5|) . This allows 
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one to define a regular finite initial value problem at space-like infinity. 

Under the evolution defined by the extended reduced system the set 1° 
evolves into a set X = ] — 1, l[xX°, which represents a boundary of the physical 
space-time. This cylinder at space-like infinity is defined solely in terms of 
conformal geometry and the general conformal field equations. 

In the given coordinates, the sets which represent near space-like infinity 
the conformal boundary at null infinity are at a finite location. They 'touch' 
the set X at certain critical sets = {±1} x which can be regarded as 
the two components of the boundary of X and, simultaneously, as boundaries of 

. Due to the peculiarities of the construction the solution is determined on 
the closure X = X U X~ U X+ ~ [—1,1] x X° of X uniquely by the data on S and 
there is no freedom to prescribe boundary data on X. 

This setting discloses the structure which decides on the asymptotic smooth- 
ness of the fields. At the critical sets X* occurs a break-down of the hyperbolicity 
of the reduced equations. As explained in section 15.31 a subtle interplay of this 
degeneracy with the structure of the initial data near X" , which is mediated by 
certain transport equations on X, turns out critical for the smoothness of the 
conformal structure at null infinity. This peculiar situation is not suggested by 
general PDE theory, it is a specific feature of Einstein's theory, the geometric 
nature of the field equations, and general properties of conformal structures. 

The transport equations, which are linear hyperbolic equations interior to 
X, allow one to calculate the coefficients of the Taylor series of the solution 
at X from data implied on X" by the Cauchy data on S (cf . 15.821) . Near X* this 
series can be interpreted as an asymptotic expansion. The coefficients are 
smooth functions on X which can be calculated order by order by following an 
algorithmic procedure. 

It turns out that the coefficients develop in general logarithmic singu- 
larities at the critical sets X*. This behaviour foreshadows a possible non- 
smoothness of the conformal structure at null infinity. In the linearized setting 
it follows that the logarithmic singularities are transported along the generators 
of null infinity ("39 ). In the non-linear case their effect on the conformal struc- 
ture at null infinity is not under control yet, however, the solutions are unlikely 
to be better behaved than in the linear case. 

The evidence obtained so far suggests cases which range from conformal 
structures of high differentiability to ones with low smoothness at null infinity. 
The non-smoothness may take the form of polylogarithmic expansions in terms 
of expressions c(l — t)'' log' (1 — r). Here r is a coordinate with r ^ 1 from 
below on J'~^ , the coefficients c = c{p, t) are smooth functions of a coordinate p 
along the null generators and suitable angular coordinates t on the set of null 
generators of and fc, j are non-negative integers. If k is small enough Sachs 
peeling fails and Penrose compactification results in weak differentiability. 

Can one 'loose physics' if one insists on extensions which are smooth at 
null infinity ? This certainly would be true if the coefficients c would encode 
important physical information. The discussion of the regularity conditions in 
section 15.51 suggests that at low orders the coefficients are determined by the 
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data in an arbitrarily small neighbourhood of space-like infinity. By the results 
on the constraints referred to above these data seem to be rather arbitrary, only 
weakly related to 'the system' characterized by the data on S, and thus of little 
relevance. As described in the following the situation is more complicated at 
higher orders, depends then in a more subtle way on the non-linearity of the 
equations, and requires further study. 

Since the coefficients can be calculated at arbitrary orders, we expect 
that this analysis will also allow us to describe in detail the relations between 
physical concepts defined at space-like infinity and concepts defined on null 
infinity (ADM resp. Bondi linear and angular momentum, etc.). The behaviour 
of the fields at the sets I* is also critical for the possibility to reduce the BMS 
group to the Poincare group (cf. 03). Thus, the precise understanding of the 
behaviour of the fields near the critical sets should provide us with a rather 
complete physical picture. 

Eventually one would like to make statements on the smoothness of the 
solution space-time at null infinity in terms of properties of the initial data. 
Thus one needs to control how the behaviour of the asymptotic expansion at 
the critical sets depends on the structure of the initial data and to derive reg- 
ularity conditions on the initial data which are necessary and sufficient for the 
smoothness of the coefficients at the critical sets. 

The information on the coefficients which is available so far has been used to 
derive conditions which are necessary for the regularity (|37|). The derivation 
of the complete condition is still difhcult because of the algebraic complexities 
of the calculations involved. 

Recently J. A. Valiente Kroon obtained with the help of an algebraic com- 
puter program complete and explicit expressions at higher order which are point- 
ing at the possibility that asymptotic staticity (or, if the time reflection sym- 
metry is dropped, asymptotic stationarity) may play a decisive role in deriving 
sufficient regularity conditions ([T^, |Z1I)- 

We are thus led to revisit the static vacuum solutions (sections 14.21 |H1 and 
[71). Because of the loss of hyperbolicity of the conformal field equations at 
I*, it is not clear whether the smoothness of the conformal structure at null 
infinity observed for static and stationary vacuum solution can be understood 
as resulting from the possible regularity of the extended solutions at the critical 
sets. We show that for static solutions our setting is smooth, in fact real analytic, 
in a neighbourhood of the set U X U X+ . 

This narrows down the range in which the final regularity condition is to be 
found. We know that asymptotic staticity is sufficient and that the conditions 
found in |37| . which are implied by asymptotic staticity, are necessary for the 
regularity of the asymptotic expansion on X. There is still the possibility that 
the final condition ' fizzles out' and depends on the specific data but we expect 
to arrive at the end at a definite, geometric condition. 

To assess how restrictive such conditions would be, it is instructive to con- 
sider the results by Chrusciel and Delay in jf6j . They allow us to conclude that 
there exist large classes of solutions to the constraints, which are essentially 
arbitrary on given compact subsets of the initial hypersurface S, whose evolu- 
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tions in time admit asymptotic expansions at with coefficients wliidi extend 
smootlily to the critical sets I* up to a given or at all orders. 

So far we ignored a question which is of central importance for gravitational 
wave astronomy: can the replacement of an asymptotically flat extension by 
another one result in a drastic change of the wave signal at null infinity ? If 
that were the case, it would be hard to see how specific physical processes could 
be identified in the wave forms calculated at J'^ . 

There should be characteristics of wave signals which are specific to 'the 
system' and which are stable under changes of the asymptotically flat extension 
if these extensions are restricted to 'reasonable' classes. This problem should be 
amenable to analytical and numerical investigations and we expect our analysis 
to contribute to its solution. In fact, it appears that with the regularity condi- 
tions mentioned above the field equations themselves hint at a first 'reasonable' 
class of asymptotically flat extensions. 

2 Conformal field equations 

Our analysis of the gravitational field near space- like and null infinity relies on a 
certain conformal representation of Einstein's vacuum field equations referred to 
as the general conformal field equations. We give a short introduction to these 
equations and point out various facts about the equations and the underlying 
mathematical structures which will be important in the following. For deriva- 
tions, detailed arguments, and further background material such as the theory 
of normal conformal Cartan connections, which is the natural home of many of 
the concepts used in the following, we refer the reader to the original article |36| 
and the survey article (38) . 

The aim is to discuss a solution (A/, g) to Einsteins vacuum field equation 

RAa] = 0, (2.1) 

in terms of a suitably chosen conformal factor Q and the conformal metric 
g — g. Denoting by V the Levi-Civita connection of g, the transformation law 
of the Ricci tensor under the conformal rescaling above takes in four dimensions 
the form 

Rup[9]^Rup[9]~^^.^pe-g,pg^' (^1 6 - ^ Va 6 6^ . (2.2) 

If O is considered here as being given, equation (|2.1|) implies with H2.2|l an 
equation for g with a similar principal part as (|2.1|l . 

As explained in the next section, we will mainly be interested in the be- 
haviour of the field in space-time domains where 0^0. Because the right 
hand side of H2.2() is formally singular in this limit, an abstract discussion of the 
solutions near the sets {Q = 0} becomes very delicate. It will be seen, however, 
that under suitable assumptions on the initial data for the field and with an 
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appropriate behaviour of the conformal factor the right hand side of H2.2|l can 
attain smooth Umits. This result is obtained by a more sophisticated use of the 
behaviour of the fields and the equations under transformations which preserve 
the conformal structure. 

2.1 The general conformal field equations 

In [2HI, has been obtained a system of equations which is regular in the 
sense that there occur no factors 9^^ on the right hand sides or factors 8 in 
the principal part of the equations. Its unknowns are 0, g, and fields derived 
from them such as the rescaled conformal Weyl tensor jki = jki- 
These have been used to derive various results about the asymptotic behaviour 
of solutions to the Einstein equations. The specific behaviour of the conformal 
fields near space-like infinity discussed in the next sections asks, however, for a 
particularly careful analysis of the equations and the gauge conditions. It turns 
out that this is considerably simplified by making use of the full freedom offered 
by conformal structures. 

A Weyl connection for the conformal structure defined by 5 is a torsion free 
connection V which satisfies 

^ p9tiy ^ Ip9t^tyi (2-3) 

with some 1-form /p. It is distinguished by the fact that it preserves the con- 
formal structure of g (and thus of g). If a frame {ek\k=o,i.2,3 is conformal at 
a given point p in the sense that it satisfies there g{ej, Ck) = r^jk with some 
A > and rjjk — diag (1, — 1, — 1, — 1), then it satisfies such a relation with a 
point dependent function A along a given curve 7 through p if it is parallely 
transported along 7 with respect to the connection V. In particular, if the 
1-form fp is closed the connection V is locally the Levi-Civita connection of a 
metric in the conformal class. 

Assuming under g — > g = 0^ 5 the transformation law 

the defining property (|2.3|) is expressed in terms of the metric g equivalently by 
^ p9pv = —'i'fpfj^iv where V denotes the Levi-Civita connection of g. It follows 
from (|2.3|) that the connection V defines with the connection V the difference 
tensor V ~ V = S{f ) given by the specific expression 

S{f)p ", = 5" ,U + 5'.f,-g,, g"' fx. (2.4) 

This, in turn, can be used to specify V in terms of V and the 1-form fp. The 
three connections are related by 

V-V = S'(/), V ~V ^ S{Q-^ dQ), V-V = 5(/). (2.5) 

Important for us will also be the 1-form 

c;^ = e/^-e/^ + v^e. (2.6) 
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The decomposition 

i?^ = 2 {g^ Lp], - .g,[A ip] ^ + C"^ .Ap, (2.7) 

of the curvature tensor of V in terms of the trace free conformal Weyl tensor 
i,\p and the Schouten tensor 

_ 1 1 
Lfiy ^ -Rfj,u - —Rg^iiy, (2.8) 

which carries the information about the Ricci tensor R^^ = Rp i^^p^, , has an 
analogue for V which takes the form 

i?^ ,Ap = 2 {g^ [A Lp]. - . L[Ap] - g,[x Lpj + C'' ,Ap, (2.9) 

where 111 

= - R(^u) ~ ^ R[piu\ - -^RQfj^i^, (2.10) 

contains the information about the Ricci tensor R^^^ — Rp ^pi, and the Ricci 
scalar R = g^^" R^^ . The tensors (|2.8ll and H2.10|l are related by 

— fp, fu + gpu fx — Lpi, — Lpy. (2.11) 



To take care of the specific direction dependence of the various fields near 
space-like infinity it is convenient to express the conformal field equations in 
terms of a suitably chosen orthonormal frame field. Let {efe}fe=o,i,2,3 be a frame 
field satisfying gik = g{ei, Cfe) = rjik, denote by Vfc and Vk the covariant deriva- 
tive with respect to V and V in the direction of , and define the connection 
coefficients Ti^ k and Ti^ k of V and V in this frame field by ViCfe = k^k 
and V^efe = T^ ^ ^ respectively. Then Ti^ k ^^i^ k + S^ i fk + S^ k h-gik 3^' /; 
with fk — fp.£^ k, where k — < dx^ , Ck > denote the frame coefhcients with 
respect to an as yet unspecified coordinate system x^, fi = 0, 1,2,3. We note 
that ft^ jti'' k because Ti ^ k gji + T^ ^ ; g^k = 0. 

If all tensor fields (except the Cfc) are expressed in terms of the frame field 
and the corresponding connection coefficients, the conformal field equations used 
in the following are written as equations for the unknown 

u=ie^k, t^k, 4-fc, W'jki), (2.12) 

and are given by 

[ep,eg] = itp\-t,'p)ei, (2.13) 

Bp{Tq j) — eq{Tp j) — Ffc ji^p q~^q 'p)~^^p k^^q j~^q fcTp j (2.14) 
= 2 {g' [p Lq]j - g' j L[pq] - g^[p L^j + 9 M^' jpq, 
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Vp Lqj — \7 q Lpj — di jpq, (2.15) 

V,PF%fei-0. (2.16) 

The square brackets in the first equation denote the commutator of vector 
fields. The connection V, which appears in the last equation, can be expressed 
by the relations given above in terms of V and fk- The last equation, referred 
to in the following as the Bianchi equation, is in a sense the core of the system. 
It is obtained from the contracted vacuum Bianchi identity C ^Ap = by 
using the specific conformal identity ri~^V^tC^^Ap = (17"^ C ^Ap)- The 
first three equations are then essentially the structural equations of the theory 
of normal conformal Cartan connections. 

No equations are given so far for the fields Q and = 8 /fc + V^G. They 
reflect the conformal gauge freedom artificially introduced here into Einstein's 
field equations. These fields cannot be prescribed quite arbitrarily. For solution 
for which the limit — > is meaningful the latter should imply dk — > Vfc0. 

The theory of normal conformal Cartan connections associates with each 
conformal structure a distinguished class of curves which provides a useful way 
of dealing with the gauge freedom. A conformal geodesic for {M,g) is a curve 
x{t) in M which solves, together with a 1-form / = /(r) along it, the system 
of ODE'S 

iVixr + S{f)x%x^xP = 0, (2.17) 

(Vi/). - i /p ^(/)a . = iA. i\ (2.18) 

where S{f) and L are given by (|2.4|) and H2.8|l with g replaced by g. For any 
given metric in the conformal class there are more conformal geodesies than 
metric geodesies because for given initial data x^ £ M, e Tx,M, g T*^M 
there exists a unique solution x{t), f{T) to (|2.17|) . H2.18|l near satisfying for 
given e M 

xin) = x,, i(n)=i„ = (2.19) 

The sign of g{x,x) is preserved near a;* but not its modulus. 

Conformal geodesies admit, unlike metric geodesies, general fractional linear 
maps as parameter transformations. They are conformal invariants. Denote by 
b a smooth 1-form field. Then, if x{t), f{T) solve the conformal geodesies 
equations H2.17|l . (|2.18() . the pair x[t), f[T) — b\x{T) solves the same equations 
with V replaced by the connection V = V -I- S{b) and L by L, i.e. the curve 
x{t), and in particular its parameter t, are independent of the Weyl connection 
in the conformal class which is used to write the equations (cf. flU]). 

Let there be given a smooth congruence of conformal geodesies which covers 
an open set U oi M such that the associated 1-forms / define a smooth field on 
U . Denote by V the torsion free connection on U which has with the connection 
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V difference tensor V — V = S{f) and denote by L the tensor H2.10|l derived 
from V. Comparing with (|2.11l) . we find tfiat equations H2.17|l . H2.18II can be 
written 

Vii = 0, Lx^x^ = 0. (2.20) 
Let efc be a frame field satisfying along the congruence 

Vi Cfc = Vi efc+ < /, Cfc > i+ < /, i > Cfe - g{x, ek)gHf, ■) = 0- (2.21) 

Suppose that 5 is a hypersurface which is transverse to the congruence, meets 
each of the curves exactly once, and on which g{ei,ek) = Q^rjik with some 
function 8* > 0. It follows that g{ei,ek) = rjik on U with a function Q 
which satisfies 

Vie = e<i,/>, ei^-e,. (2.22) 

The observations above allow us to construct a special gauge for the confor- 
mal equations. Let 5 be a space-like hypersurface in the given vacuum solution 
{M,g). We choose on 5 a positive 'conformal factor' 6*, a frame field e^*, and 
a 1-form /* such that Ol g{ei^,,ek*) — rjik and eo* is orthogonal to S. Then 
there exists through each point e 5 a unique conformal geodesic {x{t), /(t)) 
with r = on 5* which satisfies there the initial conditions x — eg*, f — f*- 

If all data are smooth these curves define in some neighbourhood U oi S 
a smooth caustic free congruence which covers U. Furthermore, / defines a 
smooth 1-form on U which supplies a Weyl connection V as described above. 
A smooth frame field Ck and the related conformal factor Q are then obtained 
on U by solving (|2.21() . H2.22|l for given initial data = e^*, 9 = 6* on S. 
The frame field is orthonormal for the metric g — g. Dragging along local 
coordinates x", a — 1,2,3, on S with the congruence and setting = r we 
obtain a coordinate system. This gauge is characterized on U by the explicit 
gauge conditions 

x^eo = dr, fo^fc^O, Lok = 0. (2.23) 

Coordinates, a frame field, and a conformal factor as above are said to define a 
conformal Gauss gauge. Since metric Gauss systems are well known to quickly 
develop caustics, it may be mentioned that conformal Gauss systems can cover 
large space-time domains in a regular fashion (cf. j4t)|). 

To obtain a closed system for all the fields entering equations (|2.13|l - (|2.16|l . 
we could now supplement the latter by equations which are implied for the 
fields 6 and dk in a conformal Gauss gauge. It turns out that such a gauge im- 
plies quite simple ordinary differential equations along the conformal geodesies 
defining the gauge, it holds in fact 

do=0, e = do, da = 0, a =1,2, 3, 

where the dot denotes differentiation with respect to the parameter r. 
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Thus, the fields Q and dk given by a conformal Gauss can be determined in 
our situation exphcitly f |36p: If g is a solution to Einstein's vacuum equations 
1^2.1]) . the fields O and dk are given by the explicit expressions 

e^e* (i + r <f,x>,+^el J (2.24) 

= e. (^1 + r < /, X >, +^ /)) J , 

do = e, da = < ea* >, a = 1,2,3, (2.25) 

where denotes the contravariant version of g and the quantities with a sub- 
script star are considered as constant along the conformal geodesies and given 
by their values on S. 

Assuming for Q and dk the expressions and (1223), equations - 

l|2.16|l provide a complete system for u. In spite of the fact that we use a special 
gauge, we refer to this system as the general conformal field equations to indicate 
that they employ the full gauge freedom preserving conformal structures. 

Equally important for us are the facts that the expression H2.24|l offers the 
possibility to control in a conformal Gauss gauge the location of the set where 
9^0 and that (|??^ . imply with the relation dfc = 6/^ + V^O in 

sufficiently regular situations that 

Vfcev'^e^o as e^o. (2.26) 



2.2 Spinor version of the general conformal field equations 

Writing the conformal equations in the spin frame formalism leads to various 
algebraic simplifications. We introduce here only the basic notions of this for- 
malism and refer the reader to |HJ for a comprehensive introduction. It should 
be noted, however, that our notation does not completely agree with that of 
j()l| . In particular, if a specific frame field is used it will always be pointed out 
in the text but not be indicated by gothic indices. 

Starting with the orthonormal frame introduced above we define null frame 
vector fields caa' = aa' Cfc with constant van der Waerden symbols cr'^ aa' 
(here and in the following indices A,B,.. ., A' , B' , . . . take values and 1 and 
the summation rule is assumed) such that 

eoo' ^ + ^3)' eii' = (eo - 63), 

1 . • ^ 1 / 

eoi' = (ei - 162), eio' = (ei + z 62). 

Then cqo', en' are real and eoi', eio' are complex (conjugate) null vector fields 
and their scalar products are given by g{eAA' , ^cc) = ^Ac ^A'C" where eACj 
^A'c : e'^*^, e'^ denote the anti-symmetric spinor fields with eoi = eo'i' = = 
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gO 1 _ 2^ rpjjg latter are also used to raise and lower spinor indices according 
to the rules = e"^^ Xb and Xb = X^ eab so that ca^ ^ f-AC e^*^ denotes 
the Kronecker spinor (similar rules hold for primed indices) . 

If connection coefficients Taa' cc are defined by writing Vg^^, ecc = 
Taa' cc sbb'j the spinor connection coefficients are given by Taa' ^ c = 
^ ^AA' CE' and one has 

BE' -r S ^ B' I f. B' ^ B 

i AA' cc — -L AA' C^C + ^ AA' C ^C ■ 

Here it is observed, as usual, that the relative order of primed and unprimed 
indices is irrelevant and that under complex conjugation primed indices are 
converted into unprimed indices and vice versa. Covariant derivatives of spinor 
fields are now given by similar rules as in the case of the standard frame for- 
malism. Writing Vaa' for Vg^^, we have e.g. for a spinor field X^ b ^ 

DD' X^ B^ = euD' {X^ B ) + ^DD' ^ F X^ b 

r F C C" vA F' 

— i DD' B ^ F + i DD' F' B 

We have similar rules for the connection V and its associated connection coef- 
ficients T AA' BC and it holds 

TcCAB = ^CCBAi ^CCAB — ^CCAB — ^AC IbC, (2.27) 

so that Tec ^ F — fee gives the 1-form relating the connection V to V. 

The general conformal field equations are now written as equations for the 
unknowns 

GAA', TaA'BC ^AA'BB', 4>abcd- (2.28) 

Here Qaa'bb' is the spinor representation of the tensor field X^j. It admits a 
decomposition of the form 

1 

QaA'BB' = ^AA'BB' ~ ^ ^AB ^A'B' + ^AB ^A'B' + ^A'B' ^AB, 

where ^aa'bb' = ^bb'AA' = ^aa'BB' represents the trace- free part of the 
tensor ^ R{jk) provided by V while R = g^'^ Rjk is the Ricci scalar and the last 
two terms, with = ^(ah)) represent the anti-symmetric tensor j R[jk]- The 
symmetric spinor field 4>abcd = 4'(abcd) represents the rescaled conformal 
Weyl tensor and is related to the latter by 

WaA'BB'CCDD' = —'t'ABCD ^A'B' ^C D' — 4>A'B'C'D' ^AB ^CD- 

The general conformal field equations in the order H2.13|l , H2.14|l , H2.15|l , H2.16|l 
now take the form 

[eBB' , ecc] = i^BB' cc - Tec bb' ) baa' , (2.29) 
ecc{TDD'^B)-eDD'{tcc''B) (2.30) 
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F t\ a , F -r A fi F'fi A 

— ^CC D^FD' B + ^ DD' C FC B " ^ CC D' ^ DF' B 

+^ DD' C'i-CF' B+^CC Fi-DD' B ~ ^ DD' F i- CC B 

= —Qbd'cc (-d^ + Obc'dd' ec^ + Qip^ bcd ^cd', 

y BB'&CC'AA' — yAA'QcC'BB' — {(j^EABC^E'C'^A' B' + (t>E' A' B'C'^EC^Ab), 

(2.31) 

A'4>ABCF^0, (2.32) 

with the fields 8, dj^A' as given above. The simple form (|2.32f) of the spinor 
version of the Bianchi equation will be useful for us. 

2.3 The reduced conformal field equations 

The conformal Gauss gauge is not only distinguished by the fact that it is 
provided by the conformal structure itself and supplies explicit information on 8 
and dfc, but also by a remarkable simplicity of the resulting evolution equations. 
Setting p = in 12.13|l - H2.16|l and observing the gauge conditions H2.23|l we 
obtain 

dre%^-tq\d'u (2.33) 
dr f / , = -f fe % f / + .9% i« + .9* , ^ go - 5,0 L,' + QW' ,0, , (2.34) 

drLqj+fq''oLk3=d^W' ,0q, (2.35) 

While the first three equations are then ordinary differential equations along 
the conformal geodesies, we still have to deduce a suitable evolution system from 
the last equation. The Bianchi equation represents an overdetermined system 
of 16 equations for the 10 independent components of the rescaled conformal 
Weyl tensor. It implies a system of wave equations for jki which could be 
used as the evolution system. For the application studied in this article it turns 
out important, however, to use the first order system. 

There are various ways of extracting from the Bianchi equations symmet- 
ric hyperbolic evolution systems but these are most easily found in the spin 
frame formalism. With the spinor field r"^'^ = eo"^ ^o' ^ + ei' ^ the gauge 
conditions (|2.23|) can be written 

r^^'eAA' = V2a„ t^^'Taa'^'c^O, t^^' 8aa'bb' - 0. (2.36) 
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Observing taa' = ^ and its complex conjugate version, one obtains 
from H2.27|l and (I2.36|) the relation t^'^' Tcc'AB = -ta fsc and thus 

^CC'AB — ^CC'AB — ^AC ^DD'EB C' , (2.37) 

which shows with H2.27|l that Tqc'AB can be expressed in our gauge in terms 
of Tcc'AB and vice versa. 

Transvecting equations (|2.29|) . (|2.30fl . H2.31|l suitably with t^^' thus gives 
the system of ODE's 

V2dr e'' cc = -Tec BB' T^^' e'' aa', (2.38) 

V2dr TdD' ^B + (TdD' C TfC ^B+ toD' ^' C f CF' ^ b) t^*^' (2.39) 
= QbC'DD- T^^' +0 0-^ BCD D', 

ecC'AA' + {TaA' B QcC'FB' + ^AA' ^' B' Occ BF')t''^' (2.40) 
= —d^^ {4'EABC^E'C''''^ A' + (pE'A'B'C'^ECTA^ )• 

We set now Aabca' = A'4>abcf- Equation (|2.32|) is then equivalent to 
— h^ABCD = A ABC A' To ^ ■ On the other hand we have the decomposition 
Aabcd = \abcd) - I e_D(c Aab)f with irreducible parts 

A{ABCD} = --^ ^ABCD — '^^'Di^D^ <t>ABC)F) ^ ^ABF ^ = 'D^^ (j)ABEF , 

(2.41) 

where P — t^^ V aa' = \/2Veo and Vab = '''{A^ ^ b)A' denote covariant 
directional derivative operators such that I?oo — ~Voi', fii — Vio', and I?oi = 
DiQ = --i= , (cf . , 1(36, for more details of the underlying space-spinor 

formalism) . 

In a Cauchy problem one will in general assume eo to be the future directed 
normal to the initial hypersurface S. The operators "Dab then involve only 
differentiation in directions tangent to S and the equations Aab f ^ — are 
interior equations on S. They represent the six real constraint equations implied 
on S by the Bianchi equation. 

For a symmetric spinor field '4'Ai...Ak we define its (independent) essential 
components by ipj — tp(Ai...Ak}j ^ where < j < k and the brackets with subscript 
j indicate that j of the indices in the brackets are set equal to 1 while the others 
are set equal to 0. 

The five equations A(abcd) = for the components of (f>ABCD contain the 
operator P. Multiplying by suitable binomial cocfhcicnts (and considering the 
frame and connection coefficients as given), wc find that the system 
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has the following properties. If (j) denotes the transpose of the C^-valued 'vector' 
(^0: 4>ii4>2i4>ii 4>4), it takes the form 

with a C^- valued function H{x,(j)) and 5 x 5-niatrix- valued functions A^^ which 
are hermitian, i.e. '^A'^ — A'-'-, and for which there exists at each point a covector 

such that A^^^fj, is positive definite. The system H2.42II is thus symmetric 
hyperbolic (03, see also and the references given there). 

While the constraints implied on a given space-like hypersurface are deter- 
mined uniquely, there is a large freedom to select useful evolution systems. In 
fact, any system of the form 

= 2aAoooi', 
= (c - d) Aooii' - 2 a Aoooo' , 

0={c + d) Aoiii' -(c-d) Aooio' , (2.43) 
= 2eAiiii- ~ {c + d) Aoiio', 
= -2eAino', 

with a, c, e > and —(2 e + c) < d < 2 a + c, is symmetric hyperbolic (the system 
(|2.42|l occurs here as a special case). We note that only the characteristics of 
these systems which are null hypersurfaces are of physical significance. 

Equations ifT^ . ifO^ . lf05|l . respectively equations (EHOI), 
combined with a choice of (|2.43|) . will be referred to in the following as the 
(general) reduced conformal field equations. Solutions to these equations solve 
in fact the complete system (|2.12|) . (|2.13|) . (|2.15() . (|2.16() if the solution admits a 
Cauchy hypersurface on which the latter equations hold ([321). In other words, 
propagation by the reduced field equations preserves the constraints. 

2.4 The conformal constraints 

To analyse solutions to the conformal field equations in the context of a Cauchy 
problem one needs to study the conformal constraints implied on a space-like 
initial hypersurface S. It will be convenient to discuss the evolution equations in 
terms of a conformal factor Q which differs on S from the one used to analyse the 
constraints. We thus assume Einstein's equations H2.1|l . a conformal rescaling 

9 = ^^9, (2.44) 

with a positive conformal factor f2, and denote again the Levi-Civita connection 
of g by V. It is also convenient to derive the conformal constraints from the 
metric conformal field equations. The latter are written in terms of the unknown 
fields g,n, S= jVf.Vn +^Rn, L^^ as in EHJ, and W^' = l^-^C^ 
and are given by equation (|2.7|) . with y\p = y\p, and the equations 

2VlS (2.45) 
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\/^\/^VL = -nL^^ + Sg^j,^, (2.46) 
which are obtained by rewriting the trace and the trace free part of (|2.2I) . 

5 -L^.V^l^, (2.47) 



VaV-Vp^a. = V^l^iy^.Ap, (2.48) 
which both can be obtained as integrabihty conditions of l|2.46|l . and 

VpW'^ = 0. (2.49) 

In these equations the Ricci scalar R is considered as the conformal gauge source 
function. Its choice, which is completely arbitrary in local studies, controls 
together with the initial data and dil on S the evolution of the conformal 
scaling. 

To derive the constraints induced by these equations on S we choose a g- 
orthonormal frame field {ek}k=o,i.2,3 near S such that n = cq is ^-normal to 
S, write = Ve^, Vfe Cj = Ffc'je;, and express all fields (except the e^) 
and equations in terms of this frame. We assume that indices a, 5, c, .. from the 
beginning of the alphabet take values 1,2,3 and that the summation convention 
also holds for these indices. The inner metric h induced by g on 5 is then given 
by hab = g{ea, Sb) = —Sab and the second fundamental form by 

Xab = gC^e^n, Cb) = Ta % gjb = -Ta ° 6- 

We set E = Vo and W*^^^ = \ W^.^apf-"^ Xp- If the tensor fields 

are projected orthogonally into S and expressed in terms of the frame {ea}a=i,2,3 
on S, they are given by (the left hand sides of) 

Lab, La = LaQ, 
Wabcd = Wabcd, Wab = WaObO, W*^ EE W*^^^, Wabc = WaObc, 

respectively and satisfy the relations 

i? = 6L^^ = 6(ioo + ^a"), (2.50) 

Wabcd ^ -2{ha[cWd]b + hb[dWc]a}, W^^ be = Wabc, W^^ = Wabc <^d'"' , 

Wab = Wba, Wa" =0, W^b = wl^, ° = 0, 

Wabc = -Wacb, w" ac = 0, W[abc] = 0, 
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where indices are moved with hab and Cabc is totally antisymmetric with £123 = 1. 
The tensors Wab and represent the n-electric and the n-magnetic part of 
W jki on S respectively. 

Equation (|2.7|) implies Gauss' and Codazzi's equation on S 

Tafc = - Wab + Lab + " hab + Xc " Xab - Xac Xb " , (2-51) 
Db Xca - Dc Xba = ^ Wabc + hab Lc - hac Lb, (2.52) 

while equations H2.45|l - H2.49|l imply the following interior equations which only 
involve derivatives in the directions of e^, a = 1, 2, 3, tangent to S 

2nS~T.'^ -Dan 0^^1 = 0, (2.53) 

DaDbn = -j:Xab-nLab + Shab, (2.54) 

Da^ = Xa''Dcn-nLa, (2.55) 
DaS = -D^nLba-^La, (2.56) 

Da Lbc - Db Lac = D^fl Wecab ~ S Wcab ~ Xac Lb + Xhc La, (2-57) 

Da Lb - DbLa = DTl Wcab + Xa "Lbc ~ Xb "Lac, (2.58) 

D^Wcab ^ X" aWbc- X" bWac, (2.59) 

D^ Wab = X''" Wabc, (2.60) 

where Vab denotes the Ricci tensor of hab- These equations can be read as 
conformal constraints for the fields 

n, S, S, 

hab, Xab, La, Lab, Wab, W^^. 

Alternatively, if a 'physical' solution hab, Xab to the vacuum constraints is given 
and a conformal factor f2 and functions S, R have been chosen, which are 
gauge dependent functions at our disposal, the equations above can be used to 
calculate 

S, L^„ 

from the conformal first and second fundamental forms hab, Xab of S, which are 
related to the physical data by 

hab = hab, Xab = n{Xab + ^hab)- (2-61) 

The equations above will be discussed in more detail in section^ 
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3 Asymptotic simplicity 



To characterize the fall-off behaviour of asymptotically flat solutions at null 
infinity in terms of geometric concepts Penrose introduced the notion of asymp- 
totic simplicity {^^, [HOIj cf. also j^Jj for further discussions and references). 

Definition 3.1 A smooth space-time {M,g) is called asymptotically simple if 

there exists a smooth, oriented, time- oriented, causal space-time {M,g) and a 

smooth function ^ on A4 such that: 

(i) M. is a manifold with boundary J , 

(a) n> on M\J andn ^0, dQ ^ on J, 

(Hi) there exists an embedding ^ of M onto ^{M) — A4\J^ which is conformal 
such that Jl^ ^^^*g — g- 

(iv) each null geodesic of(A4,g) acquires two distinct endpoints on J . 

We note that only the conformal class of iyM.,g) enters the definition and it 
is only the conformal structure of (A^, 5) which is determined here. The set J is 
referred to as the conformal boundary of {A4 , g) at null infinity. This definition 
is the mathematical basis for the 

Penrose proposal: Far fields of isolated gravitating systems behave like that of 
asymptotically simple space-times in the sense that they can be smoothly extended 
to null infinity, as indicated above, after suitable conformal rescalings. 

Since gravitational fields are governed by Einstein's equations, the proposal 
suggests a sharp characterization of the fall-off behaviour implied by the field 
equations in terms of the purely geometrical definition (|3.1|l . 

We will be interested in the following in solutions to Einstein's vacuum 
field equations (|2.1(l which satisfy the conditions of definition H3.1|l (or suitable 
generalizations). The two assumptions have important consequences for the 
structure of {M,g). We shall only quote those which will be used in the following 
discussion and refer the reader for further information to the references given 
above. 

If the vacuum field equations hold near J^, the latter defines a smooth null 
hypersurface of Ai (cf. equation (|2.45l) '). It splits into two components, J'~^ 
and , which are generated by the past and future endpoints of null geodesies 
in A4 and are thus called future and past null infinity (or scri ±) respectively. 
Each of J^^ is ruled by null generators, each set of null generators has topology 
S*^, and have the topology of M x S*^. For the applications one will have 
to relax the conditions of definition H3.1|l . In particular condition (iv), which 
is important to obtain the result about the topology of must be replaced 
by a different completeness condition if one wants to discuss space-times with 
black holes. 

One of the main difficulties in developing a well defined concept of outgoing 
radiation in the time-like cut model is related to the fact that there exists in 
general no distinguished null direction field along the time-like boundary T. In 
contrast, the null generators of J''^ define a unique causal direction field on 
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which is represented by V"^"* n. It turns out that the field (j)ABCD 
on where o'^ denotes a spinor field satisfying o"^ o"^ = — V'^'^ £7 on and 
4'ABCD the rescaled conformal Weyl spinor field, has a natural interpretation as 
the outgoing radiation field. Further important physical concepts can be asso- 
ciated with the hypersurface or subsets of it and questions of interpretation 
have been extensively analysed (cf. jj, 021 and the references given there). 

Critical however, and in fact a matter of controversy, have been the smooth- 
ness assumptions in the definition, which encode the fall-off behaviour imposed 
on the physical fields. It is far from immediate that they are in harmony with 
the fall-off behaviour imposed by the field equations. No problem arises if the 
proposal can be justified with C°° replaced by C*^ with sufficiently large integer 
k. But there is a lower threshold for the differentiability, which is not easily 
specified, at which the definition will loose much of its elegance and simplicity. 

In inO] it is assumed that 

M is of class C''+^ and g,Vl C''{M), fc > 3. (3.1) 

The conformal Weyl spinor '^abcd = ^4>abcd is then in C^^'^{Ai). Under 
the further assumption 

^ysE'^'^A'-^ABCD^Q at J+, (3.2) 

which will certainly be satisfied if fc > 4 in (|3.1|) . and the natural assumption 

the set of null generators of has topology S*^, (3.3) 

it is then shown that ^ abcd vanishes on J^'^ . The solution is thus asymptoti- 
cally fiat in the most immediate sense and the rescaled conformal Weyl spinor 
(/>ABCD extends in a continuous fashion to As a consequence, it follows 

that the space-time satisfies the Sachs peeling property (|60|. |63| . |64|') which 
says that in a suitably chosen spin frame the components of the conformal Weyl 
spinor fall-off as ^abcd = 0(f^+s+c+£'-5) (^here A, B, C, D take values 0, 1) 
along an outgoing null geodesic when its (physical) affine parameter r —^ oo at 
J+. 

Remarkable as it is that such a conclusion can be drawn for the spin-2 nature 
of the field '^abcd and its governing field equation '^abcf — 0, there 

remains the question whether the long time evolution by the field equations is 
such that assumptions H3.1|l . (|3.2|l or the conclusion drawn from them can be 
justified. 

As discussed in the introduction, we know by now that these conditions can 
be satisfied by non-trivial solutions to the vacuum field equations. What is 
not known, however, is how the solutions satisfying these conditions are to be 
characterized in terms of their Cauchy data, whether these conditions exclude 
solutions modelling important physical phenomena, and if they do, what exactly 
goes wrong. Obviously, these questions can only be answered by analysing the 
Cauchy problem for Einstein' field equations with asymptotically flat Cauchy 
data in the large with the goal to derive sharp results on the behaviour of the 
field at null infinity. 
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The results obtained so far on the existence of solutions which admit (partial) 
smooth boundaries at null infinity make it clear that the key problem here is 
the behaviour of the fields near space-like infinity. We shall not consider any 
further the results which lead to this conclusion (cf. for a discussion and 
the relevant references) but concentrate on this particular problem. 

To begin with we have a look at the asymptotic region of interest here in 
the case of Minkowski space. If the latter is given in the form {J\4 ~ = 
Vfii'dy^ dy"), the coordinate transformation <i> : a;'' = {—yxy'^)^^y^ 

renders the metric in the domain D = {y\ < 0} = {x\ < 0} in the form 
g = Tj^i, dx^ dx^ with Vl = —x\ x^ . The metric 

g^Q^g^Tj^^dx^'dx", (3.4) 

thus extends smoothly to the domain T> of points in {x\ x^ < 0} which are 
obtained as limits of sequences in V. The point = in this set then represents 
space-like infinity for Minkowski space. With this understanding it is denoted 
by i". The hypersurfaces J ^ = \x\ x^ = 0, ±x" > 0} C I? represent parts of 
future and past null infinity of Minkowski space close to space-like infinity and 
are generated by the future and past directed null geodesies of g through i". 

Consider the Cauchy hypersurface S — {y° = 0} of Minkowski space. The 
subset iS n I? is mapped by 4> onto {x" — 0, ^ 0}. Extending the latter 
to include the point x'^ = amounts to a smooth compactification 5 — *■ 5 = 
5U {i} ~ such that the point i with coordinates x^ —Q represents space-like 
infinity with respect the metric induced on iS by ^. The distinction of space-like 
infinity i with respect to a Cauchy data set and space-like infinity with respect 
to the solution space-time will become important and much clearer later on. 

Denote by hap and Xap the metric and the extrinsic curvature induced by 
^ on 5. A global representation of the conformal structure induced on S by 
haf} is obtained by using a slightly different conformal rescaling than before. 
Set h' = h with n' = 2(1 + |y|2)-i where |y| = y^WF^^Wf^W?- In 
terms of standard spherical coordinates 0, (j) on S and the coordinate x defined 
by cot ^ = |y|, < x < tt, the conformal metric takes the form h' = — (c?x^ + 
sin^xdcr^) of the standard metric on the 3-sphere and extends smoothly to 
the point i, which is given by the coordinate value x = and distinguished 
by the property that = 0, df2 = 0, Hess^l = ch! , with c ^ at «. Here 
d(7^ = dO^ + sin^ 9 d(jP' denotes the standard line element on the 2-sphere . 

Since Xq/3 = and we are free to choose S = in H2.61|l . we get x'ap = 0- 
By the formulas given in the previous section one can derive from the con- 
formal Minkowski data {S,h',x') and a suitable choice of initial data for the 
gauge defining fields H2.24II . 12.25|l a conformal initial data set for the reduced 
conformal field equations. These allow us then to recover the well known confor- 
mal embedding of Minkowski space into the Einstein cosmos ( 60 ) as a smooth 
solution to the regular conformal field equations ([3HI, |40| ) . 

We would like to control what happens if the conformal Minkowski data are 
subject to finite perturbations. Under which assumptions will the correspond- 
ing solutions preserve asymptotic simplicity ? This or the apparently simpler 
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question under which conditions the solutions will preserve near space-like in- 
finity a reasonable amount of smoothness of the conformal boundary cannot be 
answered by immediate applications of the conformal field equations. The rea- 
son is that the conformal data will not be smooth at the point i. The structure 
of the conformal initial data as well as the initial value problem for the confor- 
mal field equations near space-like infinity thus require a careful and detailed 
analysis. This will be carried out to some extent in the next section. 

4 Asymptotically flat Cauchy data 

As indicated in the case of Minkowski space above we will assume that the data 
for the conformal field equations are given on a 3-dimensional manifold S = 
iSU {i} which is obtained from a 'physical' 3-manifold S with an asymptotically 
flat end by adjoining a point i which represents space-like infinity. The data 
hab, Xab On S are thought as being obtained from the physical data hab, Xab 
by equations (|2.61() with suitable choices of and E such that all fields extend 
with an appropriate behaviour (to be specified more precisely below) to i and 
E(i) ^ 0, n > on S, = 0, Da^ii) = 0, DaDb^ii) = -2 hab, where we 
assume the notation of subsection l2.4l 

The constraint equations H2.51|l - H2.60|l contain analogues of the vacuum 
constraints. The form of these equations suggests a solution procedure which 
does not require us to go back to the physical data. By taking the trace of 
equation (|2.51|) . using (|2.53|) and the trace of H2.54|l . and writing A;i = Da I?", 
one gets 

n^r^-ArtAhil + G Dan D'^n - 4 r! E xc " + f^'((Xc ')' - Xac x""), (4.1) 

where r denotes the Ricci scalar of h. With 9 = this equation takes the 
form of Lichnerowicz' equation 

{Ah^\r)e^-U {{Xc'f~XabX'''') + \o^^Xc''- (4.2) 
By taking the trace of (|2.52l) and using 12.55|l one gets 

D\n-'' Xbc) = X& ' - 2 (4,3) 

Equations (|4.2I) and H4.3|l correspond to the Hamiltonian and the momentum 
constraint respectively. Assuming now 

Xa ° = and (the choice of gauge) E = on 5, (4-4) 

which imply Xa " = 0, equations 14. 2() and H4.3|) suggest to proceed as follows: 
(i) prescribe h on S and solve the equation 13" tpab = for a symmetric /i-trace 
free tensor field ipab on S, (ii) solve equation 14.211 with Xab = S^'^^'ab for a 
positive function 9, i.e. solve 

{Ah-\r)e = \9-'xabX''\ e>0. (4.5) 
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The fields = ^, hab, and Xab — ^"^^Pab then solve H4.1|l and H4.3|l . If it is 
required that 

pe^l, Vafc = 0(^) asp^O, (4.6) 
P 

where p{p) denotes near i the /i-distance of a point p from i, the fields /lab and 
Xab related by (j2.61(l to hab and Xab satisfy the vacuum constraints and are 
asymptotically flat (IHll)- 

Using the conformal constraints to determine the remaining conformal fields 
one gets 

s^^A^n + ^nir + xab x"") = ^ D,n D-n, (4.7) 

La = ^D-nXca, (4.8) 
Lab -^L,^hab = -^(^DaDtn-^ A/^fl hab^ , (4.9) 

Loo^lR~L,^ = ^R-^{r + Xab x"") , (4.10) 



Wq6 = (Oa DbQ - ^ Ahfl habj - ^ (xac Xb " - ^ Xce X""" Kb + Sab 



(4.11) 



<b = -^D,Xe(,a^b)'\ (4.12) 

where we set Sab = fab ~ hab- The differential identities ifT^ - which 
are not needed to get these expressions, will be then also be satisfied (cf. ISTjl. 

In view of conditions H4.6|l most of these fields will in general be singular at 
i. One will have Wab = 0{r~^) near i whenever the ADM mass m of the initial 
data set hab-, Xab does not vanish. Controlling the time evolution of these data 
requires a careful analysis of these singularities. As a simplifying hypothesis 
we assume, as in |37| . that the data are time reflection symmetric and define a 
smooth conformal structure, i.e. 

hab e C°^{S), Xab = 0. (4.13) 

We note that much of the following discussion can be extended to more general 
data such as those considered in [^S] and the more general class of data discussed 
in [23, which includes the stationary data. 

The Ricci scalar R is at our disposal. With i? = | r one gets on S 

Kb = -^ (^Da Dbfl-^ hab^ hab, (4.14) 

Loo = 0, Loa = 0, wlb^O, (4.15) 
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Wab = - ^ - - Ai.n hab + ^ Sabj ■ (4.16) 

In spite of this simplification the crucial problem is still present; one finds that 
Wab = 0{p~^) near i if m 7^ (cf. 1(2211)) ■ 

4.1 Time reflection symmetric asymptotically flat Cauchy 
data. 

To allow for more flexibility in the following analysis, we also want to admit 
non-trivial cases with vanishing or negative mass. The positive mass theorem 
(pb]) then tells us that we must allow for non-compact S. This will create 
no problems because we are interested only in the behaviour of the fields near 
space-like infinity. 

Let x", a = 1,2,3, denote /i-normal coordinates defined on some convex 
open normal neighbourhood U oi i so that with h = hab{x'^) dx°' dx^ 

x''{i) = Q, x" habix") = -x" 5ab on U. (4.17) 

All equ ations of this subsection will be written in these coordinates. We set 

\x\ = \/ Sab x°- and T = \x\'^ = Sab x"" x^ so that 

/i''''i:>aTi:)bT = -4T, (4.18) 

and 

T(z) = 0, DaT{i) = Q, DaD,T{l) = -2hac. (4.19) 

By taking derivatives of H4.18|l and using (|4.19|l one obtains 

DaDbD,Tit)^0, DaDbD,DdT{i)^~^raicd)b[h]{i), (4.20) 
where the curvature tensor of h is given by 

rabcd[h] = 2{ha[cld]b + hb[dlc]a} 

with lab[h] = Vabih] — J r[h] hab because dim(5) = 3. Proceeding further in this 
way on can determine an expansion of T in terms curvature terms at i. The 
relations above imply in particular 

(A,,T + 6)(i) = 0, DaiAhT + 6)ii)^0, Da DbiAhT + 6){i) = ^rab{i). 

(4.21) 

Equation (|4.5fl and the first of equations (|4.6(l can be combined under our 
assumptions into 

{Ah-lr)e = AnS,, 
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where in the coordinates x"" the symbol denotes the Dirac-measure with weight 
1 at = 0. In a neighbourhood of i there exists then a representation Q = 
j^ + W with functions U, W which satisfy 

{^^^-\r) =A^5,, {Ah-\r)W = Q near z, (4.22) 



and 

= VK(z) = -, (4.23) 

where m denotes the ADM-mass of the solution. The functions C/, W are 
analytic onU \i h is analytic ( ) and smooth if h is ( ) . 

The function tr = ^ is characterized uniquely by the conditions that it is 
smooth, satisfies the equation (A?i — \ r) a~^l'^ = At: and the relations 

= 0, Daa{i) = 0, Da Db <j{t) = -2 hab, (4.24) 

hold, which follow from H4.19|l and (|4.23() . If a' is another function satisfying 
these conditions, then cr' = Tt/'-^ with U' ^ 1 + 0{\x\) by and U' S 

C°°{U) by the results of The function / = ct^^/^ - ct'-^/^ ^^^^^ g^iygg 

(A,, - i r) / = and it follows that / G C°°(Z^) and \x\ f = U - U' € C^{U). 
Expanding / and U — U' in terms of spherical harmonics it follows from the 
last equation that / vanishes at i at any order. Since / satisfies the conformal 
Laplace equation it follows that / = on by theorem 17.2.6 of [HJ- This 
implies that a' — a owhi. 

The first of equations (|4.22|l can be rewritten in the form 

2 T DaU + (Aft T + 6) [/ - 2 T (A,, - i r) C/ = 0. (4.25) 

8 

This allows us to determine from (|4.24|) recursively an asymptotic expansion 
of t7, which is convergent if h is real analytic. The Hadamard parametrix 
construction is based on an ansatz 

oo 

U = ^Up TP, (4.26) 

by which the calculation of U is reduced to an ODE problem. The functions Up 
are obtained recursively by 



[/o = expJ i £\a„T + 6)^ 



(4p-2)T^ Jo Uo 
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where the integration is performed in terms of the affine parameter p = T 2 = |a;| 
along the geodesies emanating from i. It follows that 

[/(«) = !, DaU{l)=Q, DaDbU{l) = ^lab[h], (4-27) 

which implies 

DaDbD,a{i) = Q, DaDbDcDd<7{i)^2{h,dlab + hablcd)- (4.28) 

Given h and the solution W of the conformal Laplace equation in (|4.22|) , the 
considerations above show us how to determine an expansion of the function 

n = 0-^=(^ + w] - 77^^^-FT7T7, (4.29) 



/a 



in terms of p at all orders. Corresponding expansions can be obtained for the 
conformal data ||1T|) . (gUHJl, (HUnj). 

While U is thus seen to be determined locally by the metric h, the function 
W carries non-local information. Cases where 9"a W{i) — for all multiindices 
a = {a^,a^,a^) S with |a| = + + < N for some non-negative 
integer A'^ or for N = oo will also be of interest in the following. In the latter 
case we have in fact (\51\) 

W = near i. (4.30) 
For convenience this case will be referred to as the massless case. 
A rescaling 

with a smooth positive factor satisfying i9(i) — 1, leaves h = il^^ h unchanged 
but implies changes 



\x'\ 



where \x'\ is defined in terms of /I'-normal coordinates x"' as described above. 
Due to the conformal covariance of the operator on the left hand sides of equa- 
tions (|4.22l) . relations (|4.22l) . H4.23|l will then also hold with all fields replaced 
by the primed fields. 

To reduce this freedom it has been assumed in |37j that the metric h is given 
near i on 5 in the cn-gauge. By definition, this conformal gauge is satisfied by 
h if there exists a 1-form at i such that the following holds. If x(r), 1{t) solve 
the conformal geodesic equations (with respect to h) with a;(0) = i, = 
and h{x,x) = 1 at i, then a frame Cq which is /i-orthonormal at i and satisfies 
Dx Ga = (with D — D = S{1)), stays /i-orthonormal near i. This gauge can be 
achieved without restrictions on the mass and fixes the scaling uniquely up to 
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a positive real number and a 1-form given at i. It admits an easy discussion of 
limits where m — > 0. 

If ni > 0, it is convenient to set above -d = — W. It follows then that 
W = f , whence = (A^/ - | r[h']) W = -§ r[h']. Thus, if m > 0, we can 
always assume h to be given such that 

r[h]=0, n = j^--^—^ with VM=y. (4.31) 

In this gauge the function a satisfies near i the equation (cr^^/^) = inSi, 
which implies by H4.24|l 

2as = DaaD"-a with s=-Aha, (4.32) 

3 

(note that an analogous equation holds with a replaced by fi). Equation H4.32|l 
implies in turn together with H4.24|l the Poisson equation above. 

For later reference we note the form of the conformal Schwarzschild data in 
this gauge. In isotropic coordinates the Schwarzschild line element is given by 

dS^ = ( ^'"^^^[ X dt^ -{l + m/2rf {df^+PUa^). 

Expressing the initial data /i, x induced on {t = 0} in terms of the coordinate 
p = l/r, one finds that x = and h = fl^^ h with 

h = ~idp' + p'da'), 1] = ^_^-^, (4.33) 
so that a = T — resp. U = 1. The metric h also satisfies a cn-gauge. 

4.2 Static asymptotically flat Cauchy data. 

Static solutions to the vacuum field equations can be written in the form 

g = v'^ df + h, 

with t-independent negative definite metric h and t-independent norm v = 
y/g{K, K) > of the time-like Killing field K ~ dt- With the g-unit normal of 
a slice {t = const.} being given by n = v^^ K and the associated orthogonal 
projector by /i^ = g^j,'^ — fif^ny , one gets for the second fundamental form on 
this slice 

because it is symmetric by h being hypersurface orthogonal while the second 
term is anti-symmetric by the Killing equation. The solutions are thus time 
reflection symmetric. 
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For these solutions the vacuum field equations are equivalent to the require- 
ment that the static vacuum field equations 



rab[h] ^^DaDb V, A,- v ^ 0, (4.34) 

hold on one and thus on any slice {t ~ const.}. In harmonic coordinates these 
equations become elliptic and h and v thus real analytic. 

We consider solutions h, v to equations (|4.34|l with non- vanishing ADM-mass 
which are given on a 3-manifold S which is mapped by suitable coordinates i° 
difFeomorphically to R'^ \ K, where K is a closed ball in M^. We assume that 
h satisfies in these coordinates the usual condition of asymptotic flatness and 
V 1 as |i| ^ oo. The work in ^ (cf. also (54) for a strengthening of this 
result) then implies that the conformal structure defined by h can be extended 
analytically to space-like infinity. The physical 3-metric h therefore belongs to 
the class of data considered above. 

For such solutions it follows from the discussion in j7_ that the gauge H4.31|l 
is achieved if any of the equivalent equations 




2 1-v 



m 1 



(4.35) 



holds. The set S — S U {i} can then be endowed with a differential structure 
such that the metric h — il^ h extends as a real analytic metric to i. We 
shall consider in the following /i-normal coordinates as in (|4.17|l such that the 
functions a{x'^), habix"^) are then real analytic on U. The first of the static 
vacuum field equations H4.34|l then implies 

= T.ab = Da Db cr - s hab + cr (l - fi a) rab[h], (4.36) 

where s is defined as in H4.32|l . The second of equations (|4.34|) implies r[h] = 
and can thus be read as a conformally covariant Laplace equation for v. Using 
the transformation rule for this equation and observing (|4.35|) . we find that it 
transforms into 

O^{Ah-lr[h]){ev)^{Ah-lr[h]){0^m) on S, 
o o 

and is thus satisfied by our assumption r[h] = 0. 

We shall repeat some of the considerations of [SJ in the present conformal 
gauge. The fact that solutions to the conformal static field equations are real 
analytic and can be extended by analyticity into the complex domain allows us 
to use some very concise arguments. We note that the statements obtained here 
can also be obtained by recursive arguments. This will become important if some 
of the following considerations are to be transferred to C°° or C'^ situations. 

From 14.36|l one gets 

Dc T,ab^ DcDaDbCT- DcShab + <j{l~ ^J. <j) -De ^afc + (1 - 2 ^ Cr) (T Tab- 
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With the Bianchi identity, which takes in the present gauge the form D° rab = 0, 
follow the integrability conditions 

0=^D'^,a^DaS+{l-fi a) Tab D'> a, (4.37) 

and 

G = D^,^a]b + \D-'i:a[cha]b (4.38) 

= cr { (1 - ^ cr) D^c,ra\b " M (2 D^^a ra]b + D'^a r^ic ha]b)} ■ 

Equation H4.36|l thus implies an expression for the Cotton tensor, which is given 
in the present gauge by b^ca = £'[c''a]&j and for its dualized version, which is 
given by 

bab = l bacd eb = — ^ [D, a r^a tb"' ~\ r^e a t^a ■ (4.39) 
2 1 — /I cr 2 

It follows that 

I?a(2 S - i?c rr) = (7 D^S,a - 2 E,a, 

which shows that equation (|4.32l) is a consequence of equations (|4.24|l and H4.36|l 
and that the latter contain the complete information of the conformal static field 
equations. 

Let Ca = e'^a^ajc, a — 1,2,3, now denote the /i-orthonormal frame field on 
U which is parallely transported along the /i-geodesics through i and satisfies 
e°a = (J'^a at i. In the following we assume all tensor fields, except the frame 
field Ba and the coframe field a'^ dual to it, to be expressed in term of this frame 
field and set £)a = i'ea. The coefficients of h are then given by /lab = — ^ab- 
Any analytic tensor field ra/.'.'.ap' on Y has an expansion of the form (cf. |37| 'l 

Ta'';.;.a':' (^) = E ^ • ■ • (^<=^ • • • Ta^/.::a':')(*), 

k>a 

(where the summation rule ignores whether indices are bold face or not). 

We want to discuss how expansions of this form can be obtained for the 
fields 

cr, S, Tab, 

which are provided by the solutions to the conformal static field equations. 
Once these fields are known, the coefficients of the I-forms = a'^bdx^, 
which provide the coordinate expression of the metric by the relation h = 
—5ac<y^b<y'^ddx^dx'^, and the connection coefficients Fa'^c with respect to 
ea can be obtained from the structural equations in polar coordinates (cf. |5()j ) 

(pcT^bipx)) =6''b + pTc^dipx) x'^ a" bipx), 

dp 
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-r- ipTa'' eipx) a'^ b{px)) = pr" eda{px)x'^ a"^ h{px). 
a p 

For this purpose we consider the data 

Ca,...aibc - n{D^^ ... Dai '^bc)(j), (4.40) 

where TZ. means 'trace free symmetric part of. These data have the following 
interpretation. Since solutions to the conformal static field equations are real 
analytic in the given coordinates x"' , all the fields considered above can be 
extended into a complex domain U' C C'^ which comprises U as the subset of 
real points. The subset M — {T — 0} of W , where we denote by T again the 
analytic extension of the real function denoted before by the same symbol, then 
defines the cone which is generated by the complex null geodesies C D O 9 
Q ^ x°'{C) = C,x1 € W through i, where a;° 7^ is constant with habX^x^ = 
at i. On M the field D'^T dx" = ~2x°- dx<^ is tangent to the null generators of 
M ■ The derivatives of Tab x^ with respect to C at i are given by the complex 
numbers 

X^ . . . X^ X^X^Dq^^ . . . Dai^hci.^^ 

= i^-i^" . . . ^H^H^ . . . l^Da, . . . Da,b^ rcDEFii) 

where the term on the left hand side is rewritten on the right hand side in space 
spinor notation and it is used that x^^ = ax'^ — l"^ with some spinor 
because a;" is a null vector. Allowing x'^ to vary over the null cone at i, i.e. 
allowing to vary over P^(C), we can extract from the numbers above the real 
quantities 

CApBp...AiBiCDEF = ^'(Ap Sp • • -DaiBi rcDEF){i), (4-41) 

which are equivalent to (|4.4()|) . Giving the data (|4.4()|l is thus equivalent to 
giving rab(C x^) x^ where varies over a cut of the complex null cone at i or 
to giving, up to a scaling, the restriction of Vab D'^T D^T to A/". The data H4.40|l 
are in one-to-one correspondence to the multipole moments considered in 

We consider now the Bianchi identity D'^ Tab = and equation (|4.38() . In 
space spinor notation they combine into the concise form 

{l^pa)DA "^TBCDE = 2 p te^bcd Da) ""fr. (4.42) 

Note that the contraction and symmetrization on the right hand side project out 
precisely the information contained in Tab D^T D^T while the contraction which 
occurs on the left hand side prevents us from using the equation to calculate any 
of the information in H4.41|l . We use equations H4.32(l , 14.37|l in frame notation. 
By taking formal derivatives of these equations one can determine from H4.24|l 
and the data H4.4()|l all derivatives of cr, s, and Tab at «■ The complete set of data 
H4.4U|I resp. (|4.41l) is required for this and these data determine the expansion 
uniquely. This procedure has been formalized in the theory of 'exact sets of 
fields' discussed in [HJj, where equations of the type (|4.42() are considered. 
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The formulation given above suggests proving a Cauchy-Kowalevska type 
results for equations (|4.32() . (|4.37l) . (|4.42() with data prescribed on Af. Although 
the existence of the vertex at i may create some difficulties in the present case, 
this problem has much in common with the characteristic initial value problem 
for Einstein's field equations for which the existence of analytic solutions has 
been shown ([30J. At present, no decay estimates for the Cap...aibc as p ^ oo 
are available which would ensure the convergence of these series. To simplify 
the following discussion we shall assume that the series considered above do 
converge. 

We return to the coordinate formalism and show that this procedure provides 
a solution to the original equation H4.36|l . i.e. the quantity Saf, defined from the 
fields (T, s, and rab by the procedure above does vanish. We show first that Ttab = 
on Af. Since cr = on Af, this amounts to showing that mat = Da Db u — s hab 
vanishes on M. Differentiating twice the equation Dacr D°'a — 2 cr s = 0, which 
has been solved as part of the procedure above, observing that D'^ Y^cd — and 
restricting the resulting equation to Af gives the linear ODE 

D^aDc-rriab = -DaD^a rricb, 

along the null generators of Af. Observing that Da Db<J = —2 hab + 0(T), this 
ODE can be written along the null geodesies x°'{C,) = C^* considered above in 
the form 

with a smooth function Aa ~ ^a(C)- This implies the desired result. In view of 
H4.37|l . (|4.38|) it shows that we solved the problem 

D''T.ca = 0, D[cYa]b^O near i, Sab = on Af. 

The first two equations combine in space spinor notation into Da ^ ^bcde = 
with symmetric spinor field "Sabcd. Following again the arguments of jfilj . we 
conclude that T,ab = 0. 
Equation (|4.39|) implies 

D^aD^ahab^Q on V. (4.43) 

A rescaling h h' ^ 'd'^h with a positive (analytic) conformal factor gives 
a a' = d^a and bab b'^^ = -ff^^bab, whence 

D'^a D'^abab {D^'a D'ababY ^ 

^-^D^a D^a 6afc + 4 cr ^-'^D'^a D^d bab + 4 cr^ ^-^D^d D^d bab- 

This shows that H4.43|l is not conformally invariant, but it also shows that the 
relation 

D-TD''TbabW = Q. (4.44) 
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implied by H4.43|l . is conformally invariant. Using again the argument which 
ahowed us to get the quantities H4.41|l . we can translate this onto the equivalent 
relations 

niDa^---Da,hc{^)) = 0, 0,1,2,..., (4.45) 

which take in space spinor notation the form (cf. (|5.89|) ') 

DiA,B,---DA,B,bcDEF){i) = 0, p = 0, 1, 2, . . . . (4.46) 

We note that for given integer p.^, > the string of such conditions with < p < 
is conformally invariant. 

Since these condition have a particular bearing on the smoothness of grav- 
itational fields at null infinity (|31], [23) and it is not clear whether static 
equations are of a greater significance in this context than expected so far, we 
take a closer look at (|4.43|l . If we apply the operators DaDf, and DaD^Dc to 
H4.43|l and restrict the resulting equation to i, we get the relations hab{i) = 
and D(^abbc){'i') = respectively, which agree with 1)4. 45|) at the corresponding 
orders because D°-bab — 0. However, if we proceed similarly with DaDbDcDd, 
we get 

Since (|4.45|) with p = 2 can be written in the form 

D{aDbhcd){i) = y /l(a6 A/i5c<i)(i), 

the relation (|4.47|l implies in particular that Ahbcd{i) = 0. It appears that 
in general this equation cannot be deduced in the present gauge from known 
general identities and (|4.44|) alone. There will be similar such conditions at 
higher orders. While the particular form of them may depend on the conformal 
gauge, the existence of properties which go beyond (|4.45|l does not. In any case, 
these observations show that there is a gap between h satisfying the regularity 
conditions (|4.45() and h being conformally static. 

This situation is also illustrated by the following observation. If the data 
provided by h are conformally fiat in a neighbourhood of i they trivially satisfy 
conditions (|4.45|l . Without further assumptions the solution 9 to the Lichnerow- 
icz equation which relates h to the induced vacuum data h = 9^ h can still be 
quite general. However, if h is static the function 6 must be very special. 

Lemma 4.1 An asymptotically flat, static initial data set for the vacuum field 
equations with conformal metric h and positive ADM mass m is locally con- 
formally flat if and only if it satisfies near i in the gauge the equation 
Tabih] — and thus in the normal coordinates j^. _??[ ) 

-Sabdx^dx^ C/ = l, 6'=-n + ^- 

\x\ 2 

Remark: This tells us that the only asymptotically flat, static vacuum data 
which are locally conformally flat near space-like infinity are the Schwarzschild 
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data (|4.33|l . The result of (I7T|), which suggests that conformal flatness of the 
data h near i and the smoothness requirement on the functions at imply 
that the solution be asymptotically Schwarzschild, can thus be reformulated as 
saying that for the given data the smoothness requirement implies the solution 
to be asymptotically static at space-like infinity. 

Proof: By H4.38|l the solution is locally conformally flat if and only if 2 D^^a r^jjf, 
hb[cfa]dD'^'^- Applying to this equation and observing H4.36|l and again 
D[cT'a]b = 0; gcts after a contraction 

D'^a Da Tab = -3 S Tafc + CT (1 - /i ff) {hab Ted t"'^ - 3 Tac n (4.48) 

This equation can be read as an ODE along the integral curves of the vector 
field D'' cr. It follows from KTR that u"" = (2 ct a is a. unit vector field 

(with direction dependent limits at i). Because of 

2 ^ \ 



u''DaT = -\ — \ + 2U-^ D^'U DaT) <Q 



its integral curves run into i and cover in fact a (possibly small) neighbourhood 
W oil. Equation can be rewritten in the form 



''Da{T^'^rbc)=At{T^'^ rrfe), 



u 

with the matrix valued function 

Atl = ^I== {s + 2 + D'^UDa T) /^^ /i% 

+ ^^^^{hbar^^-3r\h^al 
V2cr|s| 

which is continuous on U' . This implies that Vab = on U' . The remaining 
statements follow immediately from (14.17(1 and H4.31|l . 

Remark: We note that these data may be obtained in a different form if 
locally conformally flat data are given in the cn-gauge and one asks under which 
conditions they are conformally static. The data are then of the form 

Kb = -Sab, = 9= -^ + W, AhW = 0, W{t) = ^>0. 

\x\ 2 

By a rescaling h ^ h, 9 ^-'^ 9 ^ + f ^i^^i i} = they are 

transformed into the present gauge. Assuming that these data satisfy the 
conformal static field equations and expressing the resulting equation again 
in terms of hab = —Sab one finds that the solution is static if and only if 
2W DaDbW-6 DaW DbW+2habDaW D^W ^ 0. Since Vl^ > the equation 
can be rewritten in terms oi w — W^"^ , which gives 

2wDaDbW = hab Dc w D" w. (4.49) 
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Applying £>c, multiplying with w, and using twice (|4.49|) again, we conclude that 
Da Dt, DcW — 0, whence w — k + ka x"" + kab a;° with some coefficients fc > 0, 
ka, kab- This function satisfies H4.49|l if kab — habkck'^/Ak. With ja = ka/2k 
and m — 2/^/k this gives 



2 v/l + 2ja X'^ +Jaj''XbX'>' 

with constant 

That these data are equivalent to the ones considered above is seen by rescal- 
ing with '& = W . By this one achieves W = y- For the metric ^'^ h to acquire 
the flat standard form one needs to perform a coordinate transformation which 
is given by a special conformal transformation a; — > {I o Tc o I) (x) where / de- 
notes the inversion a;" x°'{5bc x^ x'^)~^ and Tc a translation x"" ^ x°' + c° with 
suitably chosen constant c°. 

5 A regular finite initial value problem at space- 
like infinity 

In the conformal extension of Minkowski space described in section|31neighbour- 
hoods of space-like infinity, which are swept out by future complete outgoing 
and past complete incoming null geodesies, are squeezed into arbitrarily small 
neighbourhoods of the point i'^. From the point of view of the causal structure 
it is natural to indicate space- like infinity by a point. The discussion in section 
01 shows, however, that in general i*^ cannot be a regular point of any smooth 
conformal extension. The condition for an extension to to be C°° (under our 
assumption H4.13|l ') is that the data are massless in the sense of 14.31)11 and that 
the free datum h satisfies the conditions 1)5. 89|) with = oo ( 34 , 37 ). Thus, 
smoothness at «° excludes the physically interesting cases. 

A direct discussion of the initial value problem for the conformal field equa- 
tions with initial data on an initial hypersurface S = SU {i} such that jki = 
0{p~^) at i as discussed in section 0] faces considerable technical problems. Not 
only the functional analytical treatement of a corresponding PDE problem poses 
enormous difficulties but already the choice of gauge becomes very subtle. 

The setting described below has been arrived at by attempts to describe 
the structure of the singularity as clearly as possible and to deduce from the 
conformal field equations a formulation of the PDE problem which still preserves 
'some sort of hyperbolicity' at space-like infinity. It is based on conformally 
invariant concepts so that possible singularities should be identifiable as defects 
of the conformal structure. 

In a conformal Gauss gauge based on a Cauchy hypersurface S it turns 
out that after blowing up the point i into a sphere I'^ and choosing the gauge 
suitably, one arrives at a formulation of the initial value problem near space- 
like infinity in which the data can be smoothly extended to and across . In 
that gauge also the evolution equations admit a smooth extension to space-like 
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infinity. The evolution and extension process then generates from the set X° 
a cylindrical piece of space-time boundary diffeomorphic to ] — 1, l[xX°, which 
is denoted by X. It represents space-like infinity and can be considered as a 
blow-up of the point . This boundary is neither postulated nor attached 'by 
hand'. 

In this gauge the hypersurfaces J'^ representing null infinity near space-like 
infinity are given by finite values of the coordinates which are explicitly known 
(it has to be shown, of course, that the evolution extends far enough). These 
hypersurfaces touch the cylinder T at sets diffeomorphic to 1°, which can 
be thought of as boundaries of T and of respectively. The structure of the 
conformal field equations near the critical sets appears to be the key to the 
question of asymptotic smoothness. 

It may appear odd to sqeeze space-time regions of infinite extend into ar- 
bitrarily small neighbourhoods of a point i^ and then perform a complicated 
blow-up to resolve the singularity on the initial hypersurface which has been 
created by the first step. The point of the construction is that the finiteness 
of the sets allow us to disclose, to an extent that we can put our hands on 
it, a subtle feature of the field equations which otherwise would be hidden at 
infinity (in the standard vacuum representation) or in the singularity at i^ (in 
the standard conformal rescaling). 

In the following the setting indicated above and its various implications will 
be discussed in detail. While we shall add more recent results we shall follow to 
a large extent the original article [23 ■ For derivations and details we refer the 
reader to this or the articles quoted below. 

5.1 The gauge on the initial shce and the blow-up at i 

The non-smoothness of the conformal data H4.14|l , 14.7|l , (|4.15|l , H4.16|l , H4.29|l 

at i arises from the presence of various factors p in the explicit expressions. To 
properly take care of the specific radial and angular behaviour of the various 
fields it is natural to choose the frame field in the general conformal field equa- 
tions such that the spatial vector fields Ca, a = 1, 2, 3, are tangent to the initial 
hypersurface S and one of them, 63 say, is radial. Since there is no preferred 
direction at i, this only makes sense if the frame is chosen on S such that it has 
direction dependent limits at i. This singular situation finds a well organized 
description in terms of a smooth submanifold of the bundle of frames. To dis- 
cuss the field equations in the spin frame formalism, we will consider in fact a 
submanifold Ce of the bundle of normalized spin frames over S near i. While 
the use of spinors leads to various simplifications, it should be mentioned that 
the construction could be carried out similarly in the standard frame formalism 
(cf. i38,). 

5.1.1 The construction of Ce 

Consider now 5 as a space-like Cauchy hypersurface of a 4-dimensional solution 
space-time {M, g) with induced metric h on S. Denote by SL{S) the set of spin 



38 



frames S — {Sa}a=o,i on S which are normaHzed with respect to the alternating 
form e, such that 

e((5A,fe) = CAB, £01 = 1- (5.1) 

The group 

SLi2, C) ^{t^se GL{2, C) | SAct^ st'^ d ^ esi?}, 

acts on SL{S) hy 5 5 ■ t ^ {5a b}s=o,i- The vector field r = V^cq, with 
eo the future directed unit normal of S, defines a subbundle SU{S) of SL{S) 
which is given by the spin frames in SL{S) with 

g{T, SaSa') = eA °eA' °' + ^e^' ^' = taa'- (5.2) 
It has structure group 

SU{2) - {t^ B e 5i(2, C) I TAA't^ B t^' B- = TBB'}- 

In any frame in SU{S) the vector r is given by . In the following we use 
the space spinor formalism in the notation of |36| . Using the van der Waerden 
symbols for space spinors 

^ AB „ (A ^B)A' c ^ A'c ^ 1 O Q 

CTa —CTa A'T , O AB ~ T(^g a A)A', C=L,Z,6, 

which satisfy 

h ^ ^ AB ^ B ^ B' 1^ ^BB' ^ a F ^EB' B 

iT-ab — CTaABCTb , ^A ^A' —-^TAA'T + CT APT A'T VaE , 

where 

hab ASO''' CD = ~^A{C^D)B = hABCD with hab — —Sab, 

the covering map onto the connected component 50(3) of the rotation group is 
given by 

SU{2) 3t^B^t\^a-ABt^ct'' D'Jb^'' e 50(3). 

The induced isomorphism of Lie algebras will be denoted hy ^>^. 

The covering morphism of SU{S) onto the bundle 0+(5) of positively ori- 
ented orthonormal frames on S maps the frame 6 G SU{S) onto the frame with 
vectors Cq — ea{S) — (Ja^^ Satb^ 5b' such that h{ea,eb) — hab- We use this 
map to pull back to SU{S) the /i-Levi-Civita connection form on 0+{S). Com- 
bining this with the map ^'^7^, the connection is represented by an su(2)-valued 
connection form lj^ b on SU{S). Similarly, pulling back the R^-valued solder 
form on 0+ (S) and contracting with the van der Waerden symbols results in a 
1-form a"^^ on SU {S) which is referred to as solder form on SU (5) . 

Let H denote the real horizontal vector field on SU (5) satisfying < , H > 
= eo ei or, equivalently, 

Ts{tt) H{6) - (5(0 Ti) ""'Sb' = ^{So So' - ^i'), S e SU{S). (5.3) 
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It follows that Tstin) H{St) = Ts{tt) H{5) if and only if 

teU{l)^{teSU{i)\t=(^ J,^ ),0eM}. 

The field H will essentially correspond to the 'radial' vector field mentioned 
above. 

We consider again the normal coordinates satisfying H4.17|l near i, set Be = 
{p E U \ \x{p)\ < e} with e > chosen such that the closure of Be in S is 
contained in U, and denote by {SU{Be),'n') the restriction of {SU{S),Tr) to 
Be. Let S* be in the fiber n-^ii) C SU{Be) over i. The map SU{2) 3 t 
6{t) = 6* ■ t E TT^^{i) defines a smooth parametrization of 7r^^(i). We denote 
hy ] — e,e[3 p ^ 5{p,t) £ SU{Be) the integral curve of the vector field \/2H 
satisfying S{0,t) = 5{t) and set Ce = {S{p,t) € SU{Be) \ \p\ < e, t £ SU{2)}. 
This set defines a smooth submanifold of SU{Be) which is diffeomorphic to 
] — e, e[y.SU{2). The restriction of tt to this set will be denoted by tt'. 

The symbol p, which has been introduced already in section r4.1l is used here 
for the following reason. The integral curves of \/2H through 7r~^(i) project 
onto geodesies through i with h-unit tangent vector. Thus, the projection tt' 
maps Ce onto Be- The action of ?7(1) on SU{Be) induces an action on Ce- While 

= 7T^^{i) = {p = 0} is diffeomorphic to SU{2), the fiber tt ~^{p) C Ce over a 
point p in the pimctured disk Be = Be \ {i} coincides with an orbit of U{1) in 
SU{Be) on which p — \x{p)\ and another one on which p — — \x{p)\- 

The map tt' factorises as Ce C'e Be with = Ce/U{l) diffeomorphic 
to ] — e,e[xS^- For with < < e the subsets {p = p*} of Ce are 
diffeomorphic to SU{2) and the restrictions of the map tti to these sets define 
Hopf fibrations of the form 

SU {2) 3 V2(t'' AB^^at^ 1 e CM.^- (5.4) 

The set 1^2^ (Be) (resp. n ^^{Be)) consists of two components Cg* (resp. Cf) 
on which ztp > respectively. Each of the sets CJ^ is mapped by 7r2 diffeonior- 
phically onto the punctured disk. If Be is now identified via 7r2 with Cg^ the 
manifold Be is embedded into C'e such that it acquires the set 7ri(I°) = 7r^"^(«) 
as a boundary. The set Be = BelJ 7r^^(i) — [0, e[xS'^ is a smooth manifold with 
boundary. Viewing Be again as the subset of 5 = 5\ {«}, we get an extension S 
of S which can be thought of as being obtained from S by blowing up the point 
i into a sphere. This is our desired extension of the physical initial manifold 
and the following discussion could be carried out in terms of the 3-dimensional 
manifold Be- 
lt turns out more convenient, however, to use the 4-dimensional U{1) bun- 
dle C+ = C+ UI" ^ {S € Ce \ piS) > 0} ~ [0,e[x5f/(2). It is a manifold 
with boundary smoothly embedded into SU{Be), from which it inherits various 
structures. The set Ce is conveniently parametrized by p and the paralleliz- 
able group SU{2)- The solder and the connection form on SU{Be) pull back to 
smooth 1-forms on Ce - We denote the latter again by cr"^ and Co"' b respectively. 
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Any smooth spinor field ^ on Be defines on Ce a smooth 'spinor valued function' 
which is given at 6 G Ce by the components of ^ in the frame defined by S and 
denoted (in the case of a covariant field) by £,Ai...Ak.A'i---A'j- We shall refer to 
this function as to the 'lift' of 

The structure equations induce on Ce the equations 

J AB ~A . EB ~ B . AE /c c;\ 



dti^B = -w^B A (i-^s + f7^B, (5.6) 

where 

6^ — ^ „CD . KF 

" B — 2^ BCDFCr A (7 

denotes the curvature form determined by the curvature spinor tabcdef- It 
holds 

"TABCDEF ~ I 2 '^'^SCB ~ Y2 ^BCE I ^DF + I 2 ^^^DF " '^^BDF I ^CE 

(5.7) 

where sabcd = S(abcd) is the trace free part of the Ricci tensor of h and r its 
Ricci scalar. The curvature tensor of h is given by 

TAGBHCDEF — —TABCDEF^GH ~ TgHCDEF^AB- 

and the Bianchi identity reads QD^^ sabcd — DcDf- 

We use t G SU{2, C) and = p a,s 'coordinates' on Ce- The vector field H 
tangent to Ce then takes the form \/2H = dp. Consider now the basis 

"1 K ^ ) ' ^ K 1 0^ ) ' = K -0 ' ^^'^^ 

of the Lie algebra su{2). Here is the generator of the group U (1). We denote 
by Zui, i — 0,1, 2, the KiUing vector fields generated on SU{Be) by Ui and the 
action of SU{2). These fields are tangent to 1°. We set there 

X+ = —{Zu2 + iZui), = —{Zu2 — iZui), X = —2 i Z^^, 

and extend these fields smoothly to Ce by requiring 

[H,X]^Q, [H,X±]^Q. (5.9) 

The vector fields H, X, X^, X- constitute a frame field on Ce which satisfies 
besides (|5.9|) the commutation relations 

[X, X+] ^2X+, [X, X_] = -2X_, [X+, X_] = -X. (5.10) 

The vector field iX is tangent to the fibers defined by tti . The complex vector 
fields Xj^, X^ are complex conjugates of each other such that X_ / = X^ f for 
any real- valued function /. 
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These vector fields are related to the 1-forms above by 

<(T^^,fl'> = eo(^ei^), <CT^^,X> = on Ce, (5.11) 



(5.12) 

<a;^s,F> = 0, <a;^B,X> = eo^es°-ei^es' on Ce, (5.13) 



<a;^B,X+>=eo^eB'+0(p'), < b, X_ > = -ei ^ es + O(p2), (5.14) 
as p — > 0. 

To transfer the tensor calculus on Be to Cg we define vector fields cab = C(ab) 
on Ce \ 1° by requiring 

<cr^^,ccr. > = e(c'^er.)^, ccn = cVr,ap + c+cT>^+ + c- cd^-. (5.15) 

The first condition implies that ^^(Tr') cab = ^(aTs) ^'^b' for 5 e Ce while 

the second removes the freedom for the vector fields to pick up an arbitrary 
component in the direction of X. It follows that 

c'-ab=xab, ab = - zab + ab, c~ ab = -Uab + c~ ab, (5.16) 
P P 

with smooth functions which satisfy 

c"ab = 0(p), c"oi = 0, a=l,+,-, (5.17) 

and 

XAB = \/2e(A ° es) \ VAB = -^CA^CB^, zab = -^(^A°eB°- (5.18) 
The connection coefficients with respect to cab satisfy 

7CD '^B=<U)^B,CCD>=^ 1*CD B + ICD ^ B (5.19) 

with 

7* ABCD = ^{^ACXbD + £BDXac), 701CD = 0, 7ABCD = 0{p). 

The smoothness of the 1-forms and the vector fields H, X^, X- implies that 
the vector fields pccD and the functions 

CD, pC^ CD, pC~ CD, PlCDAB, 
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extend smoothly to all of Ce- 

A smooth function F on an open subset of Ce is said to have spin weight s 

if 

X{F)^2sF (5.20) 

on this set with 2s an integer. Any spinor valued function induced by a spinor 
field on Be has a well defined spin weight, it holds e.g. 

X(j)ABCD^2{2~A-B~C-D)cj)ABCD- (5.21) 

It follows from the construction of Ce that this is also true for the functions 
considered above, it turns out that 

Xc' AB ^2(1- A- B)c' AB, X C± AS = 2 (1 - (±1) -A-B)C^AB, 

XjABCD^2{2~A-B-C-D)^ABCD for A,B,C,D^ 0,1. 

By our construction, equation 1)5. 3|l . and the formula for ea{S) given above 
the vectors ^^(Tr') {^/2H{S)) = €3(6) are tangent to and the frame ea{S{p,t)) 
is parallely propagated along the geodesies [— e,e[9 p n'{S{p,t)) through i. 
Thus we have constructed the type of frame field asked for in the beginning. 
Working on Ce has the advantage that p and \/2 H define smooth fields and the 
smoothness of the various fields considered above can easily be discussed. 

The transition from Be to Cg respectively to C^ amounts to a new choice of 
differential structure at space-like infinity. This change is reflected in the drop of 
rank of the map tt' at the set X'^ . It follows from (|5.11|) , (|5.12|) , that at points over 
i the vectors X, X± project onto the zero vector while at points in tt ^^{p) the 
real and imaginary parts of H, X^, X^ have non- vanishing projections which 
span the tangent space Tp Be if p G Be- The relations H5.11|l . H5.12|l . H5.13|l . H5.14|l 
show that the behaviour of the map tt' near X° is encoded in the behaviour of 
the solder and the connection form. 

With the structures given above we can perform tensor calculations defined 
on Be now also on Ce \ and they follow the 'usual' rules of the spin frame 
formalism. If F denotes the lift of a smooth function / on Be, the covariant 
differential Df is represented on Ce \ 1° by the invariant function DabI = 
cab{F)- In the following we shall use the same symbol for a function and its 
lift. If p,AB is the invariant function induced by a spatial spinor field p on Be 
its covariant differential is given on Ce \ 1'^ by the expression 

DabP-AD — Cab{PAd) — lAB ^ C t^ED ~ lAB ^ D PCE- 

Analogous formulas hold for covariant differentials of spinor fields of higher 
valence. 

In terms of p and t = [t^ b) on Ce and the normal coordinates x° satisfying 
H4.17|l on Be, the projection tt' has the local expression 

tt' ■- {p, t) ^ x%p, t)^pV2a''cDt^ot^i- (5.22) 
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This can be used to pull back the functions fi, U, and W, which are related 
by H4.29() . to functions of spin weight zero on Ce- The metric in (|4.13(l is built 
into our formalism and the second fundamental form lifts to a symmetric spinor 
valued function xabcd which vanishes everwhere. Using the fields 

AB, 1CDAB-, Sabcd, r, (5.23) 

given by H5.16|l . (|5.19|) . and (|5.7|) . one can determine 

DabDcd^, (5.24) 

on C+ and thus also the derived data H4.14|l . (|4.15() . H4.16|l . 

In particular, a detailed expression for the rescaled conformal Weyl spinor 
4>ABCD is obtained on C+ by using H4.16|l and (|4.29() . It takes the form 

(pABCD = (t^ABCD + 4'ABCD^ (5-25) 

where 

4>'abcd = o'"^ {D(^ab Dcd) o- + cr sabcd] (5.26) 
= ^ {U^D^^abDcd) {p^)-8pUD^ABpDcD)U} 

{2UD^abDcd)U-6D^abUDcd)U-U^ sabcd } 
is derived from a ~ p^ U^"^ and thus from the local geometry near z, while 

'I^ABCD^^ {^QUWD(abpDcd)P + UWD^abDcd){p^)] (5.27) 

+ -i {W D(ABPDCD)U -iU D(ABPDCD)W) 

-^{U D^AB DcD) W + W D(^AB DcD) U-Q D(^abU Dcd) W-UW sabcd) 

-2WD(^AB DcD) W + 6 D^abW Dcd) W + sabcd. 

is the part of the rescaled conformal Weyl spinor which depends on the non-local 
information in W and which vanishes in the massless case. Observing that 

DabP = xab, Dab Dcd{p^) = --^PliAB^ cxd)e ^ 0{p^), DabU = 0{p), 

one finds that 

^'ABCD = Oi^), <f>'lBCD = ~^-^^ABCD + 0{\), (5.28) 

where we set c'abcd = Ha ^b'' ^c^ ^d) for j = 0, . . . , 4. 
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5.1.2 Normal expansions at X" and the functions T„i^ k 

To analyse in detail the behaviour of the various fields near space-like infinity it 
is convenient to study a particular type of expansion. It will be discussed here 
for an unprimed spinor field, similar expansions hold for other fields. In terms 
of the normal coordinates x°' on Be define the radial vector field V ~ x'^ dx" ■ 
Let S* = (5* (a;") be the smooth spin frame field on Be which satisfies Dy 6* — 
on Be and coincides with the spin frame at i chosen as the starting point for 
our construction of Cg. Denote by the orthonormal frame associated with 
S* and write V = V^^ e\g. 

Suppose ^ is a smooth spinor field on Be which is given in terms of the spin 
frame field 5* by ^ = ^ i^"")- Then its Taylor expansion at i is of the 
form 

p=oo 
p=0 ^' 

(5.29) 

To determine the lift ^Ai...Ai of this field to C+ one has to observe its transfor- 
mation behaviour t^^ Ai ■ ■ ■ t^' Ai under changes of the frame 
and the fact that the pull back of the functions V"^^ are given in view of H5.22|l 

by 

l^^^(x''(p,i)) = ^/2p^(^o^''^. (5.30) 

If the expansion coefficients DspCp ■ ■ ■ DbiCi^Ai then decomposed 

into products of eab's and symmetric spinors at i, the essential components 
= (.{ai...a,)j, < j < I, < j < I, which are of spin weight 5 = 5 - j, are 
obtained as expansion of the form 

00 

where 

2p+l m 

S,j,P — ^j,P;m,k Tm " in^+j (5.32) 

TO=max{|/-2j|,i-2p} k=0 

with complex coefficients ^j.p-m.k and functions ^ k oi t as discussed below. 

We refer to this type of expansion as to the normal expansion of ^ at l'^. In 
the case considered above the lift of ^ to C+ has smooth limits at . Corre- 
sponding expansions in terms oi , k E can also be obtained for fields such 
as (j)ABCD on Ce which are given by algebraic expressions of regular fields but 
which become singular at X". 

The functions Tm-' k, arise (apart from some normalizing factors) naturally 
by the procedure indicated above. They are matrix elements of unitary repre- 
sentations 

SU{2) 3 t ^ T„,(t) =. [Tm' kit)) e SU{m + 1), 



45 



which are given by 

To = 1, ^ ,{t) = {^^ " {^^ " t^-^^ , 

j,k — . . . ,m, m = 1, 2, 3, . . . . 

The brackets with lower index now indicate synimetrization and taking 'essential 
components'. The expansions obtained above make sense under quite general 
assumptions; the functions \/ to + 1 ^ ^^[t) form a complete orthonormal set 
in the Hilbert space L^{n, SU{2)) where /i denotes the normalized Haar measure 
on SU{2). 

Using the identification of X'^ with SU (2) built into our construction, we 
consider the Tm ■' k as functions on X° and extend them as p-independent func- 
tions to Ce- The vector fields X±, X then act as left invariant vector fields and 
it holds 

XT„%- = (to-2j)T„'=,, (5.33) 



X+T.m j — Pm.jTm j-l: X^T„i j — —f3„i,j + lT.m j + 1 (5.34) 

for < fc, j < to, TO = 0, 1, 2, . . ., with P^nj — \/ i j + 1)- It follows that 
functions / with spin weight s have expansions of the form 

m 

/= ^ ^ ^ ^ Sm,k T„i ^ m.-g, (5.35) 
m>|2s| k=Q 

where the to's are even if s is an integer and odd if s is a half-integer. All 
functions considered in the following have integer spin weight. 



5.2 The regularizing gauge for the evolution equations 

To obtain definite expressions for the expansions of the data at i and because 
the terms of lower order are then simplified, it has been assumed in 37^ that the 
metric h is given in a cn-gauge near i. This will be assumed also here, though 
the discussion of the static case given below will show that this is not necessary 
for our construction. The coordinates p, t and the frame field constructed above 
depend on the choice of scaling of the metric h on S. Most important is the 
fact that = 0{p^) near it affects the definition of p in an essential way. 

In analysing the evolution of our data in time it turns out convenient to use 
a different conformal factor 8 which is related to the conformal factor Vl by 

Q = K-^VL on Ce, (5.36) 

with a function 

K^pn' with K'eC°°(Ce), k'>0, Xk! = Q, k'|io = 1. (5.37) 
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The value of k! on X° is chosen for convenience here, nothing is gained in the 
following by requiring a different (positive) boundary value for it. 

The change of the conformal factor implies a map S : (5 — > k2 (5 which maps 
the set C+ bijectively onto a smooth submanifold C* of the bundle of conformal 
spin frames over Be- We use the diffeomorphism S to carry the coordinates p 
and t and the vector fields 9p,X, to C*. The projection of C* onto Be 

will be denoted again by tt'. 

Assuming a conformal Gauss system for the evolution in time as described 
in section 12.11 the evolution of the spin frames constituting C* defines in the 
the bundle of conformal frames over the space-time manifold M. a smoothly 
embedded 5-dimensional manifold M which is again a t/(l) bundle over the 
space-time and whose projection onto M. we denote again by tt'. The manifold 
C* represents a smooth hypersurface of M . 

By pushing forward the coordinates p, t and the vector fields dp, X, X± with 
the flow of the conformal geodesies ruling Af, these structures can be extended 
to Af such that iX generates the kernel of tt'. The parameter a;° = t of the 
conformal geodesies defines a further independent coordinate with = t — 
on C* , so that the tangent vector field of this congruence can be denoted by dr- 

The reduced field equations (jTi^ . (jOi^ . lfnn|l . ifO^ (the latter special- 
ization of H2.43|l is chosen here for only definiteness) are now interpreted as 
equations on J\f by assuming that the caa' are vector fields on J\f which are 
defined at a spin frame S £ J\f hy the requirement that they project onto the 
frame defined by 5 on A^, i.e. Tg n'{eAA') = Sa Sa' , and whose X-component is 
fixed by requiring an expansion of the form 

BAA' = -^TAA' dr - A' SAB, (5.38) 

with 'spatial vectors' 

eAB = e° AB dr + ab dp + e+ ab X+ + ab X^. (5.39) 

The unknowns in the reduced field equations arc then interpreted as spinor 
valued functions on M. It can be shown that spin weights are preserved under 
the evolution by the reduced system. 

We have to express the initial data for the conformal field equations in terms 
of the new scaling. With k, the fields (|5.23|) . (|5.24|) . and the associated covariant 
derivatives (carried over to C*, observing that the local expression of S in the 
given coordinates is the identity) one gets for the curvature fields 

^3 

4>ABCD = ^ {D(^ab Dcd)^ + ^ sabcd) , (5.40) 

Qaa'CC'^'K^ Q^D(^abDcd)^+ j^rhABCD^ T^A'T^c- (5.41) 
For the frame (|5.38() . one gets by (|5.16(l 

e°AB = 0, AB = pn' XAB, (5-42) 
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AB ^ ZAB + KC^ AB-, e AB^u'VAB + KC AB- (5.43) 

For the conformal factor O we get 

e = e, = K-^n= ^ ^ .„ on c*. (5.44) 

We assume that initial data for the 1-form /, which will be related in the end 
to / by the relation / = / — satisfy 

< /, > = 0, pull back of / to C* = kT^ dn. (5.45) 

It follows then that from H2.24I) that O takes the form 

e = (^1 - ^) on A^, (5.46) 
with a function which is given by 

, '^^ ^=^ p{U + pW) iU"^ + 2pUx^^DabU ~ p'^ D^^U DabU 

y/\DjlD^\ ^ 

(5.47) 

+2 p"^ U x"^^ DabW ~ 2 p^ D^^U DabW - p'^ D^^W DabW}~^ on C* . 

Here the second member is given in the notation of section l4. II while the term on 
the right hand side is given in the notation of section lS.l.ll In (|5.4()l) and in the 
following formulas the subscripts * are saying that the corresponding functions 
are constant along the conformal geodesies. 

For dAA' we get by (|2.25|1 the explicit expression 

dAA' = TAA' Q-T^ A' dAB On A/", (5.48) 

where the dot denotes the derivative with respect to t and 

, „ f UxAB-pDABU-p^DABW \ 

= [ ^TTTW ) . ' ^'-''^ 

where the notation of section 15. 1.11 is used on the right hand side. 

If one uses (|5.42ll and (|5.43|) to write for a given smooth function p onC* 

PAB = k"^ (e^ AB dp + e+ AB X+ + ab X^) p, 

one gets with the 1-form (|5.45|l and the spatial connection coefficients (|5.19|l the 
space-time connection coefficients in the form 

^AA'CD = [ 2 ^ (''-^'^ '^'bD+ <^'BD I^'aC ) - Pl^' lABCD + -j^ ^AB I^CD j A' ■ 



(5.50) 
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Note that the Taa'bc hi the reduced equations can be expressed by (|2.37|) in 
terms of the Taa'BC- 

Most important for us is the observation that the functions given by (|5.4Uf) , 
H5.41|l . H5.42|l . (|5.43|) . H5.46|l . (|5.48|) . H5.50|) /law smooth limits as p ^ and can 
in fact be smoothly extended into the coordinate range p < 0. For the unknowns 
in the new scahng we thus obtain normal expansion in terms of non-negative 
powers of p. In particular, one has 

(1>ABCD = H^i(l>'ABCD + ^t^ABCo)^ (5-51) 

with H5.26|l . (|5.27|) on the right hand side. Pushing the expansion (|5.28|) a bit 
further and using (|5.36() one gets in the cn-gauge (in which hai, = —5ab + 0{p^)) 

(t'ABCD = -K^ 6me\g(jD (5.52) 
-p 12 {X+ Wi ABCD + 3 VKi abcd - ^- Wi abcd) 

t G) fG) - ¥ ^ 

+o(p=^). 

It is assumed here that W is an arbitrary solution to (A/j — i r) = on Bf.- 
Its normal expansion takes in the cn-gauge the form 

2 2 2p 

W = Y,pPW^ + 0{p^) = ^ / (£ W-p;2p,fc ^ p) + 0{p') 

with ^ 

m /2\ ^ 

W^0;0,0 = W^(0 = y, W^l;2.fe = (J D^,,)^W*{i), 
W2.A,k=Q ' i?(afci?c.).W'*(*). 

In the case where k' is constant the right hand side of 15.52|l provides the terms 
of a normal expansion up to the quadrupole term. If k depends on p and t the 
terms given above need to be expanded further to obtain the normal expansion. 

The transition H5.36|l to the conformal factor O corresponds to a transition 
h h' = h of the metric on Be (assuming that k' arise as a lift of a smooth 
positive function on Be with K'{i) = 1) in the sense that then il"^ h = h = 
G^^ h' . The coordinate p is then not adapted to the geometry defined by the 
metric h' . To illustrate the situation assume that h is flat. Then 

h' = -k'-^ p-^ (dp^ + p^da^) = {dr^ + da^), (5.53) 

with r = — log/9 near i. With respect to the new coordinate r, which is adapted 
to the geometry of h' , the point i is shifted to infinity but the surface measure of 
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any sphere around i remains finite and positively bounded from below. This be- 
haviour is reflected by the fact that the frame coefficient ab in H5.42|l vanishes 
while the frame coefhcients ab in H5.43|) have finite and non-vanishing limits 
at 1°. We shall keep the coordinate p because it ensures the finite coordinate 
representation of the boundary as well as the smoothness of the data near 

With the gauge defined above the functions 8 and dAA' in equations (|2.39() , 
H2.40|l are given by (|5.46|l and l|5.48|l and the finite regular initial value prob- 
lem near space-like infinity for the reduced field equations (|2.38l) . H2.39|l . (|2.4U|) . 
H2.42|l is completely determined. We write this system schematically as sys- 
tem of equations for the unknown u = {w,4>) with (j) = {4>abcd) and w = 
{eAA',i^AA'BC,QAA'BB') 01, alternatively, w = (caa' ,Taa'bc,Qaa'bb')- It 
takes the form 

drw ^ F{x,w,4>), A>'d^^.4> = H{w)4>, (5.54) 

where the ^-dependence in the first equation comes in here via the functions Q 
and dAA'- 

Important for the following is that with any choice of k satisfying jj.^'Tl j 
the functions Q and dAA' take smooth limits as p — > and can be extended 
smoothly into a range where p < 0. With the smooth extensibility of the initial 
data observed before, we find that the initial value problem for the reduced field 
equations with the data prescribed above can be extended smoothly into a range 
where p < so that the reduced equations form still a symmetric hyperbolic 
system. It may be noted finally that the congruence of conformal geodesies 
(considered as point sets) underlying our gauge does not depend on the choice 
of K, whereas the parameter t depends on it in an essential way. 

5.3 Specific properties of the regular finite initial value 
problem at space-like infinity 

The nature of the initial value problem formulated above is conveniently dis- 
cussed by considering certain subsets of R x M x SU (2) which are defined by the 
range admitted for the coordinates {T,p,t). We define 5-dimensional subsets 

= {|r| < -, < p < e, te SU{2)}, 

K 

^={\t\ < -, 0<p<e, teSU{2)}, 

K 

where — is a function of p and t. It holds then 

M =AfUj- UJ+ UJUJ- UJ+, 
with 4-dimensional submanifolds 
J'^ = {t^ ±-, < p < e, te SU{2)}, I = {\t\ < 1, p = 0, te SU{2)}, 



50 



and 3-dimcnsional submanifolds 

I± = {|t| ± 1, p = 0, te SU{2)}, X" = {r = 0, p = 0, te SU{2)}, 

where it has been observed that — ^ 1 as p ^ 0. We note that 

e>0on7V', e^O, dQ^O on J-UJ+UJ, 6 0, = on I± . 

The set C* = {r = 0, < p < e, t g SU{2)} defines a hypersurface of Af. Its 
closure in Af is given by 

C = {r 0, < p < e, te SU{2)} = C* U X°. 

Factoring out the group U{1) implies projections (denoted again by tt') onto 
subsets M X M X S"^ which are of one dimension lower than the sets above. In 
particular, J\f projects onto a set A4 which represents the 'physical space-time'. 
For convenience we will usually work with the manifolds above and use for 
them the same words as for the projections, so that Af will be referred to as the 
'physical space-time' etc. 

For suitable e > consider a smooth extension of the data given on C* 
to the set Cext — {t ~ 0, — e < p < e, t S SU{2)} and an extension of the 
functions 6, dAA' to the domain Afext = {\t\ < ^, — e < p < e, t e SU{2)}, 
so that the reduced conformal field equations (jO^ . ((2301), still 

represent a symmetric hyperbolic system of the form (|5.54|l . Then there exists a 
neighbourhood V of Cext in Afext on which there exists a unique smooth solution 
BAA', ^AA'BC (resp. Taa'bc), Qaa'bb', 4>a bcd t o our extended initial value 
problem which satisfies the gauge conditions (|2.36|) . 

It turns out, that the restriction of this solution to the set VOAf is uniquely 
determined by the data on C* . The data on C* have a unique smooth extension 
to C and it follow s from (|E3HIl, EHHll, fO^ - and (jOi^ that cc ^ as 
p —^ 0. Equations (|2.38() imply in particular 

V2dr el cc = -Tec BB' r^^' aa'- (5.55) 

It follows that cc = on V n X and as a consequence that the matrices A^^ 
in (|5.54|l are such that 

={] on X, (5.56) 

if the solution extends far enough. One can apply to the system (|5.54|l on 
subsets of V n A/" the standard method of deriving energy estimates. Without 
further information on the system the partial integration would yield contri- 
butions from boundary integrals over parts of V n X. Because of (|5.56f) these 
boundary integrals vanish and one obtains energy estimates which allow one to 
show the asserted uniqueness property. The extension above has been consid- 
ered to simplify the argument. Alternatively, the space-time M can be thought 
of as a solution of a very specific 'maximally dissipative' initial boundary value 
problem where initial data are prescribed on C and no data are prescribed on X 
because of (|5.56f) (cf. and the existence theory in 02], 103) ■ 
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In the present gauge the set X, which is generated from X° by the extension 
and evolution process, can be considered as being obtained by performing Umits 
of conformal geodesies. It represents a boundary of the space-time Af which may 
be understood as a blow-up of the point i*^. We refer to it as the cylinder at 
space-like infinity. 

Suppose that there exists on A/" a smooth solution caa', ^aa'BC (resp. 
^AA'Bc), Qaa'bb', 4>abcd of the reduced conformal field equations (|2.38() . 

(E3ni), which satisfies the gauge conditions on N and coin- 

cides on the initial hypersurface C* = {r = 0} C A/" with the data given above. 
The projections Tn^CAA') then define a frame field on Ai for which exists a 
unique smooth metric g on ^A such that g{TiT' {eAA'),Ti:' [caa')) = ^ab ^A'B'- 
Denote by D' the domain of dependence in M. with respect to g of the set 
7r'(C*) and set D = tt ~^{D'). By the discussion above we can assume that 
the closure of D in A/" contains the set I and the solution extends smoothly to 
I. It follows from the structure of the characteristics of the reduced equations, 
that the solution is determined on D uniquely by the data on C* and it follows 
from the discussion in and the fact that the data satisfy the constraints 
that the complete set (|2.29|l . (|2.30|l . H2.31|l . H2.32|l of conformal field equations 
is satisfied on D. Since Q has spin weight zero it descends to a function on M. 
and g — 0^^ g satisfies the vacuum field equations. 

The restriction to D arises here because we only considered the data on C* . 
Observing that the latter were obtained by restricting the data given on the 
initial hypersurface S to Be it is reasonable to assume that the conformal field 
equations hold everywhere on A/'UI and g defines a solution to the vacuum field 
equations on M. 

Assume u is a solution of a (possibly non-linear) hyperbolic system of partial 
differential equations of first order on some manifold. A hypersurface of this 
manifold is then called a characteristic of that system (with respect to u), if 
the system implies for some components of u non-trivial interior differential 
equations on the hypersurface. These interior equations are called transport 
equations (cf. 22 ). 

Because of (|5.56|l the set T is then a characteristic of the extended field 
equations. It is in fact of a very special type (i.e. a total characteristic), because 
the system (|5.54() reduces on 2 to an interior symmetric hyperbolic system of 
transport equations for the complete system of unknowns. Together with the 
data on 2'^ it allows us to determine u — (f , </>) on I. 

Suppose that the solution extends in a fashion to the sets J^^ . Since 
8 = 0, d<d 7^ on J'^ the sets 7t'{J'^) form (part of) the conformal boundary 
at null infinity for the vacuum solution g. Since 4>abcd is one finds Sachs 
peeling. Of course, it will be one of our main tasks to control under which 
assumptions the solutions will extend with a certain smoothness to the sets J'^ . 

As remarked before, we can expect the decision about the smoothness of 
the solution at null infinity to be made in the area where the latter 'touches' 
space-like infinity. This location has a precise meaning in the present setting. It 
is given by the critical sets , which can be considered either as boundaries of 
or as the boundary components of I. The nature of these sets is elucidated 
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by studying conformal Minkowski space in the present setting. 

We start with the hne element given by (|5.53l) and choose k' — 1. Since 
L0_= p by GH it follows that ^ {|r| = ±1, < p < e, t e SU{2)} and 
M = {\t\ < 1, < p < e} X S-^. It will be useful to express the frames 
considered in the following in terms of the specific frame 

vo = df, vi = p9p, i;± = X±. (5.57) 

The complete solution to the conformal field equations then is given by 

e\A' = 7f { " ^) ° + + ^) ' eA' vo (5.58) 
+ {eA ° eA' °' -eA' e^' '')vi-eA° sa' v+ ~ ca' ea' °' v. } 

^AA'BC = ~\'^AA' XBC, (5.59) 
&AA'BB' - 0, (5.60) 
^ABCD = 0. (5.61) 

The conformal factor and the metric g implied by e^^, are given by 

e* = p(l-T2), g* = dr^ + 2-dTdp-^—:^dp^ ~ da^. (5.62) 



With the coordinate transformation 

1 



p(l-r2)' p(l-r2)' 



(5.63) 



one gets in fact the standard Minkowski metric 5 = ^ g* = dt^ — dr^ — dcr^ 
in spherical coordinates . The flat metric corresponding to H3.4|l is given by 
n*^g = p^g* = d{Tpf - dp^ - p^da^ with n* ^ pQ* ^ p^ ~ {t p)^. For this 
metric the curves with constant coordinates p, 9, and (j) are obviously conformal 
geodesies and because of their conformal invariance it follows that the corre- 
sponding curves for g* are conformal geodesies with parameter r. Equations 
H5.63|l can be read as their parametrized version in Minkowski space. 

The metric g* given by (|5.62|) extends smoothly across null infinity but it 
has no reasonable limit at X. Its contravariant version 

= (1 - r^) dl + 2Tdr {pdp) - {pdpf - [da^, 

does extend smoothly to T. While it drops rank in the limit, it does imply a 
smooth contravariant metric on I whose covariant version = (1 — t^)^^ dr^ — 
da"^ defines a smooth conformally flat Lorentz metric on I. The coordinate 
transformation r — sin^ shows that this metric is not complete. The Killing 
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fields of Minkowski space, which are conformal Killing fields for g*, extend 
smoothly to X such that they become tangent to I, vanish there in the case of 
the translational Killing fields, and act as non-trival conformal Killing fields for 
the metric /* in the case of infinitesimal Lorentz transformations. 

The fields (jESHl, EHHl, extend smoothly to all of M. The property 

H5.5(j|l results from the fact that the fields Cqq/, e\j^, become linear dependent 
on T. Since they do not vanish there, this degeneracy does not cause any 
difficulties in the field equations. On J+ and J~ however, the field Cqq, and e*^ 
respectively vanishes. This strong degeneracy has important consequences for 
the (extended) conformal field equations. To see this, we solve the transport 
equations on X to determine the matrices A^^ on X in the general case. 

Extending the data (|On|) . ifCTTl . ifO^ . (|051) . lt53n|l . one finds that they 
agree on X°, irrespective of the choice of k' satisfying conditions of H5.37|l . with 
the implied Minkowski data. Since the extensions of the functions & and (Iaa' 
vanish on X, the transport equations for the frame, connection, and Ricci tensor 
coefficients are independent of the choice of initial data. It follows that the 
restrictions of these coefficients to X agree with those of the Minkowski data 
given above. It follows in particular that aa' = on X. Applying formally 
the operator dp to equation H5.55|l (which is part of the reduced field equations) , 
restricting to X, and observing the data dpC^ aa'\i», one finds that dpC^ aa' = 
eA ° iA' °' — ^ ea' ^' on X. Writing 

eAA'=e'AA'Vi with i==0, 1,+,-, 

and assuming the summation rule, we find that irrespective of the free datum h 
given on S and the choice of k' the fields caa', ^aa'bc, Qaa'bb' coincide at 
lowest order with the Minkowski fields above in the sense that Qaa'BB' ~ 0{p) 
and 

e-^ AA' = e*"^ AA' + aa' , ^aa'bc — ^aa'bc +^aa'bc, (5.64) 

with 

AA' - 0(p), Taa'bc = 0{p) as p 0. (5.65) 

Assuming k = w in the general case, which by H5.47|l is consistent with H5.37|l 
if e is chosen small enough, the similarity with the Minkowski case becomes 
even closer. Then Q = f Q* with proportionality factor f = which extends 
smoothly to Af such that / 1 on X. The set J'^ is given as in the Minkowski 
case above. The discussion below shows, however, that this particular choice of 
K may not always be the most useful one. 

In the general case the first deviation from the Minkowski case is found in the 
value of the rescaled conformal Weyl spinor on X. On X° it is given by (|5.52|) . 
Restricting the Bianchi equation to X and using the coefficients determined 
above one finds that 

(I^ABCD = -Qme^BCD on X. (5.66) 

The discusssion above shows that the matrices are determined on X by 
the Minkowski data and the structure of the characteristics of the evolution 
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equations for the rescaled conformal Weyl tensor agrees on T with that of the 
equations which are obtained by Unearizing the Bianchi equation on Minkowski 
space. These (overdetermined) spin-2 equations take in the gauge above the 
form 

(1 + r) - p ap^fc + X+ Vfc+i + (2 - fc) ^fe = 0, (5.67) 

(1 - r) a,Vfc+i + pSpV'fc+i + X- V^fc + (1 - fc) ^fc+i - 0, (5.68) 

where fc = 0, 1, 2, 3 and the '(/'j denote the essential components of the hnearized 
conformal Weyl spinor. 

The most conspicuous feature of these equations is the factor (1 + r) in 
H5.67|l . which vanishes on U X^, and the factor (1 — r) in H5.68|l . which 
vanishes on U X+. On these factors arise because the coordinate r 
is constant on and these sets are characteristics for the equations. By 
choosing k! differently, this degeneracy can be removed on (cf. At 
I^, however, this degeneracy cannot be removed in the present setting. Any 
symmetric hyperbolic system extracted from these equations, like e.g. 

(1 + r) Bt t/jQ - pdpipo + X+ij^i = -2 -00, 

(4 + 2 r) 9, 7^1 - 2 p V^i + X_ V^o + 3 X+^j^ - -4 t/-!, 

6 9^ V2 + 3 + 3 X+i>i = 0, 

(4 - 2 r) 9, ^3 + 2 p 9p 03 + 3 X_i/;2 + = 4 V's, 

(1 - t) dr 04 + Pdp 00 + X+'4>i = 2 04. 

must contain such factors at least in the equations for 0o and 04. Writing 
this in the form 9^ = Hip, and writing —< dr,£, >, £,p =< 9p,C >, 
^± = < ^, X± > we find 

det(^'^ Cm) = 24 {g''- (3 + g'^'' C.) 

with 

g'- ^p 6 = (1 - r') e + ^rpCr^p-p'ep-l (e+ + 

It follows that characteristics pertaining to the quadratic terms which start 
on X, stay on X and that those starting in the physical space-time never end 
on I U I~ U X+ but always run out to J'^. Most importantly however, and 
this also holds true for the general system H5.54|l . the quadratic form g^'^ £,p,£,i, 
degenerates at X* and there is a loss of real characteristics. This follows also 
directly from 

det(yl^) = onX±. (5.69) 

It appears that this loss of hyperbolicity at the critical sets is the key to the 
smoothness problem for the conformal structure at null infinity. 
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5.4 The s-jet at space-like infinity 

The relations H5.56I) and (|5.()9|l are the dominant features of the regular finite 
initial value problem at space-like infinity. The consequences of (|5.69|l are not 
deduced by the standard textbook analysis, we have to rely on the specific prop- 
erties of our problem. It turns out that a considerable amount of information on 
the behaviour of the solution near the critical sets can be obtained by exploiting 
H5.56|l . We know already that the solution is smooth in some neighbourhood of 
C m. J\f and that u can be calculated on I by solving intrinsic equations on I. 
It will be shown in the following that a full formal expansion of u in terms of p 
can be calculated on I by solving certain transport equations. 

The following notation will be convenient in the following. For p = 0, 1, 2, . . . 
and any sufficiently smooth (possibly vector- valued) function / defined on AfUT 
we write for the restriction to I of the p-th radial derivative f. The set 
of functions f'~', . . . , fP on T will be denoted by Jf (/) and referred to as the 
jet of order p of f on 2 (and similarly with 2 replaced by 2^ .) If u = (w, (p) 
is a solution of equations H5.54|l wc refer to Jj{u) (respectively J|(w), Jji^)) 
as to the s-jet of u (resp. w, (j>) of order p and to the data Jjo(w) (respectively 
Jjo{w), Jjo(0)) onX° as to the d-jet of u (resp. w, 4>) of order p. A s-jet Jj{u) 
(respectively Jj{w), (0)) of order p will be called regular on 

(or simply regular) if the corresponding functions on 2 extend smoothly to the 
critical sets 2^ . 

An initial data set on S will be called asymptotically static of order p, where 
p G NU {oo}, if its d-jet Jjo {u) coincides with the d-jet of order p of some static 
asymptotically fiat data set defined on some neighbourhood of i in S. It will be 
seen later that asymptotic staticity (of order p) is an important feature of initial 
data sets. 

Applying the operator formally to the first of equations H5.54|l and re- 
stricting to 2, one obtains for an equation of the form 

drwP = G{T,t,w°,...wP-\wP,(l)°,...<j>P-^), p=l,2,..., (5.70) 

where the right hand side is an affine function of w^. The functions do not 
appear here, because the rescaled conformal Weyl spinor occurs in the equations 
for the frame, connection, and Ricci coefficients with the factors & and dAB, 
which vanish on 2. It follows that the s-jet Jj{w) can be determined by the 
integration of an (easily solvable) linear system of ODE's, if the s-jet J^~^{u) 
and the d-jet Jjq{w) are known. 

With the notation 15.64|l the Bianchi equation can be written 

V*-'^ A' (pBCDF = -4>A'BCD, (5-71) 

where 

4>A'BCD = (t>A'(BCD) = e'^ A' Vi{(j)BCDF) - 4:t^ A' ^ {B 4>CDF)E- (5-72) 
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Then 

{V2V*^ A' <j>BCDFr = ~V2 0P ,5CD> (5-73) 
provides equations with left hand sides given by 

(1 + r) dr cj/^ + X+ 0^+1 + (2 - J - p) = . . . , (5.74) 



(1 - r) dr 0^+1 + X_ + (1 J + p) 0^^^^ = . . . , (5.75) 
where j = 0, . . . , 3, and right hand sides given by (of. (|5.fci5|l ') 



p 



t'A'BCD 



4=0, + ,- 3 = 1 



"CDF)E- 



We note that these expressions depend on J^{w) but only on Jj^^((/)). Thus, 
given these s-jets, the s-jet Jj{(f>) can be obtained by solving a linear system of 
ode's, if J^o(0) is given. Because the system is singular at the critical sets it is 
not clear a priori that Jj(0) is regular, even if J^{w) and J^~^((/)) are regular. 

To obtain more detailed information on the solutions, it is useful to consider 
a system system of second order. From (|5.71ll follows 

y EE' ^ 9ABCD — ^\'A ^ E'<PBCDE~^^a (pE'BCD, 



which is equivalent to 



^*EE' V*^^ 4iABCD = IaBCD = ~'^^*E'(a4>^ BCD), (5-'''7) 

= gBC = ^*^'^(^A'ABC- (5.78) 
While the right hand side of 

(V^;^, V-^^^' ^PabcdY = Fabcd^ (5-79) 

depends again, similar to H5.76|l . on Jj^w) and J|^^(0), the left hand side takes 
the decoupled form 

(l-r2)92^P + 2{(p-l)T-j + 2}9.</.^ + C0P~p(p-l)0P = ... (5.80) 

where the spin weight relations X = 2 (2 — j) and the Casimir operator 
C = -\ {X+ X^ + X^ X+) + i on SU{2) have been used to arrive at this 
expression. 
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The fields (j)^ have expansions 

E C ^here 0^^,^ = C \-2+„ 
with coefiicients i/)^ ^ — 'P^ q ki''')- Since the Casimir operator satisfies 

C {T2q*' q-2+]) = 9 ((? + 1) ^2 q q-2+] , 

equation (|5.8UI) imphes for (/i^^ ODE's of the form 

D(„.„,^) 4,1^ = (1 - r^) + {;3 _ « - (a + /3 + 2) r} 9, 0^^^, (5.81) 

+n{n + a + (3 + I) = . . . 

with 

Q! = j— p — 2, /? = —j — p + 2, n = ni = p + q or n = n2=j3 — (? — 1. 

The equations above allow us to calculate recursively a formal expansion of 
the solution u = (w, </)) to (|5.54|) in a series of the form 

"=E^"'''^^ (5.82) 

on 2 (note the different meanings of the supersripts p) with coefficients = 
uP{T,t) G C°°{T). In some neighbourhood of X° in Sf this series represents in 
fact the Taylor series of smooth functions and it converges near 1° if the datum 
h is real analytic. We shall try to deduce from it information on the behaviour 
of u near the critical sets. 

5.5 Behaviour of the s-jets near the critical sets. 

Because Jj{u) is regular, the integration gives a regular s-jet Jj(w). The cal- 
culation of Jj{4>) gives (in the cn-gauge and with k! = 1) the regular solution 

(t^ABCD = -{Wi 36 (1 - t2) + (18 - 3 r^} ^bcd (5.83) 

-12 (1 - rf X+ Wi el abcd + 12 (1 + rf X_ Wi abcd- 

Thus Jxiw) will again be regular. It tuns out that Jj{(t)) will not necessarily be 
regular. The integration (cn-gauge, k! — 1) gives 

4>ABCD = 4>^ABCD + 'PaBCD + 4'ABCD: 

with 

<^(isci?)o = 0' ^fABCDh = C2(^) mWi+ C3(t) m^, = 0, 
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^UbcdH = c,{T)mX+Wi, cffX^CD), - -ci(-t) mX_W^i, 



where the Ci(r) are polynomials in t of order < 



^(ABCD) 



-4J6 . (l+T)-'(l-T)*-^'^W^2;4,fcr4^-, (5.84) 
V ^-^Z fe=0 



and 



with 



4 



, = «^ (^) ^ E \/ 2 ( J ^(SFGff). T, ^ , (5.85) 



ao(T) = 2 (1 - rf K{-t) = -a4(-T), 
ai(T) = 4 (1 - rf (1 + r) K{-t) - = -«3(-r), 

a2(r) = V6 {-^-^ - 2 (1 - t)^ (1 + rf K{r)} - -a2(-r), 
(1 + r)- 



s)(l + s)5' 

While the first two terms extend smoothly to , the third term has logarithmic 
singularities at the critical sets unless the regularity condition bABCoii) = is 
satisfied (the quadrupole term W2, which looks so innocent here, reappears in 
obstructions to smoothness at higher order [71], |72|). 

It is thus clearly important to control the behaviour of the s-jets at X* at 
all orders. Equations D(^n,a.p)'^ = are well known from the theory of Jacobi 
polynomials and they have been used in to derive a certain representation 
of the solutions in terms of polynomials built from the generalized Jacobi poly- 
nomials Pi°''^\T) (ISHl)- By the overdeterminedness of the system (|5.74|l . H5.75|l 
the problem can be reduced to the integration of the functions (^q ^j , (ji^ ^j . The 
functions (jji qj 4'2qi 4'3 q can be calculated from them algebraically. 

One finds for p > 3 and q ^ p the representation 

0S,p-(l-T)P+'(l + ^r' (5.86) 

(p+i)(p+2),^p , r ^ 



4p ^^"^P* ^'P*' Jo {I + t')p-^ {1 - t')p+ 

<p = (l+T)^+'(l-^r' (-^Ip*- (5.87) 
(p+l)(p + 2) ^„ ^ f-^ dr' 



4p ^^O'f* ^-P*' Jo (l + t')p-^1 - t')p+3 

where the substcript * indicates initial data on 2^ . 
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Denoting by y^.g the column vector formed from c/jq g, 4>^iqi one obtains for 
p > 3 and < g < p — 1 

2/p,«(t) = Xp,,(T) yp,,, + ^ Xp,,(T')-i Sp,,(T') dr'j . (5.88) 

The functions Bp q are derived from the right hand sides of H5.73|l and H5.79|l 
and can thus be calculated from J^{w) and J^^'^{cf)). The matrix- valued valued 
functions Xp^q are given by 

_ ( {1+tY~^p + t) 
~ 1^ (1 - r)f-2 [p 

y _({i + tY-^ 

X ^i^^-^^") (-l)'Q3;.,,(r) \ 2<o<«-l 

^p.. - 1^ (-1). Q3;p,,(-t) Qi;p,,(-t) ; ' ^ - - ^' 

with polynomials 



of degree rii 



The solutions to the transport equations can be calculated, order by or- 
der, explicitly. The only difficulty is the calculation of the functions Bp q = 
Bp^q[Jj{w), Jj^^{(f>)] which become more and more complicated at each step. 

The most conspicuous feature of these expressions is the occurrence of loga- 
rithmic singularities at . The latter can arise, as a consequence of the evolu- 
tion process and the structure of the data, even under the strongest smoothness 
assumptions on the conformal datum h. We will have to discuss to what extent 
the occurrence of such singularities can be related to the structure of the initial 
data and whether it can be avoided by a judicious choice of the latter. 



5.6 Regularity conditions 

Expanding the integrals in (|5.86() . H5.87|l one finds 

<^o,p ~ (1 - t)p+^ (1 + t)p-^ log(l - r) + analytic in t as t ^ 1 

and a similar behaviour for t/)^ ^ as r ^ —1, unless the initial data on 2^ satisfy 
the condition 

^0,p * ~ '^4,p * ■ 

(Note that the singularities get less severe with increasing p.) This raises the 
question whether data can be given which satisfy these conditions. By a lenghty 
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recursion argument it can be shown f|37p that for given integer > the fields 
p resulting from f5.86]) . \5.87}) extend smoothly to for 2 < p < + 2 if 
and only if the free datum h satisfies the regularity condition 

D(A,B,---DA,B,bABCD){^)^0, q ^ 0, I, 2, . . . , p, . (5.89) 

By (|4.46(l these conditions are satisfied for static data with p* = oo. This 
aUows one to construct a large class of data satisfying (|5.89f) by gluing with a 
partition of unity an asymptotically flat static end to a given time reflection 
symmetric data set and solving the Lichnerowicz equation. 

Condition (|5.89|) has been observed as a regularity condition before. In [311 
has been derived under the strong assumption that the solution be massless (cf. 
H4.30|l ') a necessary and sufficient condition on h that space- like infinity can be 
represented by a regular point i*^ in a smooth conformal space-time extension 
(so that J^^ will be smooth near space-like infinity). This condition, referred 
to as radiativity condition, implies l|5.89|l . It has been shown in [l^7\ that these 
two conditions are in fact equivalent. 

The first term in (|5.88|) is polynomial and thus regular. The second term 
is not so easy to handle. If (|5.89|) is not assumed the corresponding log-terms 
will enter the integral in a non-linear way and the solution will have at 2"^ 
polyhomogeneous expansions in terms of expressions (1 =F r)*^ log"'(l T t) with 
k,j G Nq. We shall assume therefore that (|5.89|) holds with p^ — oo. 

From the expressions above it follows that the Wronskian det{Xp g) has a 
factor (1 — T^)P~^. The regularity of the integrals in (|5.88|l thus depends on 
the precise structure of the functions Bp^qlr), which get quite complicate with 
increasing p. It has been shown in j87j and j41j that Jj{u) is regular for p < 3 
if H5.89|l is satisfied with p* < 1. 

Because the functions Bp^q are getting increasingly complicated with p, J. 
A. Valiente-Kroon studied the case where h is conformalUy flat on Be with the 
help of an algebraic computer program ([7l]). In that case condition (|5.89() is 
trivially satisfied but there still exists a large class of non-trivial data for which 
h is not conformally flat outside Be- In the conformal factor (|5.44|) one has 
U = 1 on Be but W will be a non-trivial solution to the conformally covariant 
Laplace equation with m = 2W{i) ^ 0. It turns out that J^iu) is again 
regular. For Jj{u) however, logarithmic terms are observed. They come with 
certain coefficients which depend on the data. Choosing the data such that 
these coefficients vanish, still new logarithmic terms are observed for Jjiu). 
Restricting to the axially symmetric case to keep the expressions manageable, 
new logarithmic terms crop up for p — 7 and p — 8. 

The form of the conditions obtained at these orders suggests a general for- 
mula which needs to be satisfied to excluded logarithmic terms at any given 
order p (TT). If this formula is correct, all derivatives of W must vanish at i 
if the logarithmic terms are required to vanish at all orders. As a consequence 
the solution must become asymptotically Schwarzschild at i (cf. Lemma |4.1|l . 
Since W is governed on Be by an elliptic equation with analytic coefficients it 
would follow that the solution is precisely Schwarzschild near i. 
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How seriously do we need to take the singularities at ? To answer this 
question one needs to control the evolution of the field in a full neighbourhood of 
T in M.. This has not been achieved yet. However, the analysis of the linearized 
setting, which is given by the spin-2 equations (|5.t)7|) , (|5.t)8|) on Minkowski space 
in the gauge (|5.58|) . gives some insight (i39j|). 

While the functions Bp^q vanish in that case, the singularities arising from 
H5.86|l . H5.87|l do in general survive the linearization process. The analysis then 
shows that for prescribed integer j the function 

p-i -, 

p'=0 ^ 

extends to a function of class on A4 , if one chooses p> j + 6 m the expansion 
above. Here "0^ i p' = 0, 1, . . .p — 1, are understood as p-independent functions 
on Ai UX, which agree on J with the s-jet Jj~^(V') (defined by equations (|5.()7|) . 
(ESHll)- Note that the sum above provides the first terms of an asymptotic 
expansion of the solution at . 

It follows that the solution will extend smoothly to all of Ai if the linearized 
version of (|5.89|l is satisfied with p^, = oo. If the condition is satisfied only 
with some finite > 2 but violated at p = + 1, the solution will develop a 
logarithmic singularity at X* which will be transported along the null generators 
of so that the solution will be only in CP*~^(7V(). While it remains to be 
seen whether the solutions to the non-linear equations admit similar asymptotic 
expansions at JT"^, the discussion shows clearly that the regularity of the s-jets 
Jjiu) is a prerequisite for the smooth extensibility of the solutions to J'^ . 

If the solutions to the non-linear equations show a singular behaviour on 
J'^ as indicated above, does it refer to something 'real' or to a failure of the 
gauge ? If the underlying conformal structure where smooth at null infinity, the 
conformal geodesies should pass through where 6^0 and the 1-form, the 
V-parallely transported frame, and therefore also the rescaled conformal Weyl 
spinor in that frame should be represented by smooth functions of r along the 
conformal geodesies because these as well as their natural parameter t depend 
only on the conformal structure. Singularities as indicated above therefore refer 
to intrinsic features of the underlying conformal structure. 

The results of (|7I]) show first of all that the regularity condition H5.89|l with 
p» = oo are not sufficient for the regularity of Jjiu), p = 0, 1, 2, . . .. It appears 
that the Lichnerowicz equation, which breaks the conformal invariance by fixing 
the scaling of the physcial metric h = fi"^ h, does play a role in the smoothness 
of the conformal structure at null infinity. This is remarkable because it shows 
that besides the local condition H5.89|l there are other conditions to be observed 
which are 'not so local'. However, the Lichnerowicz equation is introduced only 
as a device to reduce the problem of solving the underdetermined elliptic system 
of constraints to an elliptic problem. The results of ^B], [201, [21] exploit the 
underdetermincdness of the contraints in quite a different way. They teach us 
to be careful with the words 'local' and 'global' in the present context. 
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The main purpose of calculating Jjiu) for the first few p is to get an insight 
into H5.88|l which would allow us to control the behaviour of Jjiu) near m 
dependence of the data given on S. One may speculate that the results above are 
telling us that asymptotic staticity, or more generally asymptotic stationarity, 
at space-like infinity is of more importance in the present context than expected 
so far. Recent generalizations of the calculations in ([7J) to non-conformally 
flat data seem to support this view (|72|'). 

This raises the question whether the setting proposed in ||37j is for static 
solutions as smooth as one would expect. This is far from obvious because of 
the loss of hyperbolicity at the critical sets. Giving an answer to this question 
for general static solutions will be the purpose of the following chapters. 



6 Conformal extensions of static vacuum space- 
times 

For static asymptotically flat vacuum solutions with positive ADM mass we 
shall construct in the following a conformal extension which will include null 
infinity and will also allow us to discuss the cylinder at space-like infinity. The 
extension will be defined in terms of explicitly given coordinates and conformal 
rescaling. In section[3will it be shown that it coincides with the extension (not 
the coordinates etc.) as defined in section [^31 

Because one expects usually 'not much to happen at space-like infinity' for 
static asymptotically flat solutions, one may wonder why the detailed discussion 
of the fields near space-like infinity should be so complicated. An obvious reason 
is that a gauge which is chosen to discuss space-like and null infinity must 
introduce a 'time dependence', it cannot be adapted to a Killing field whose 
flow lines run out to time-like inflnity. However, the main reason is that the 
static field equations play an important role in discussing the regularity of the 
field near the critical sets; we will have to make extensive use of them. 

The static vacuum solution is assumed in the form 

g^v^ dt^ + h, 

with V = v{x'^), h = hab{x'^) dx"" dx^ and a conformal factor SI — ri(x'^), where we 
assume /i-normal coordinates which satisfy (|4.17l) and the conformal gauge 
which achieves (|4.31ll on the set M x Z^, where U — {\x\ < p*} with a sufficiently 
smaU p* > 0. We set 

x" 1 



x\^, e'' = f- = — T-i/2^°T for > 0, p - 



Y^ix^^y 

\ a=0 



Coordinates ^ = 2,3, on the sphere — {\x\ — 1} can be used to 
parametrize e° and we write then e° = e°'{4>^) and de" = e"" ^^d^p^. For 
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convenience the coordinates tp^ will be assumed in the following to be real 
analytic. If x"^ = pe°'{ilj^), the metric h takes the form 

h = —d + fP' k, 

with (p-dependent) 2-metrics 

k = kAcdip^ dip'-' = hac{pe'')de°' de", 

on the spheres p = const. > 0. For p — > the metric k approaches the stan- 
dard line clement da^ = —k{0,tl)^) on the 2-dimensional unit sphere in the 
coordinates tp^. 

We write now = t and a;°' = f , a;^ = p, x^' = and consider the map 
$ : a;''' ^ x^{x^') defined by 

p 

p(l-r) 

It follows that the four differentials 

dx" = ((l-f)rfp-prff)e" + p(l -f)de", (6.2) 



[{vn){pe-) {vn){p{i-f)e-))'^''^ {vn){p{i-f)e-)'^^^^' ^^-^^ 

with 

p 

are independent for < f < 1 and < p < p* and we can consider the as 
smooth coordinates on an open neighbourhood of space-like infinity in {t> 0}. 

For s > we set 

^.._ (^»)(^e-) U{se^)-sf 

To indicate the different arguments replacing s in this and other functions 
of s e" or of s and ip^, we write out the argument replacing s explicitly but 
suppress the dependence on e" or tp^. Thus h{s) will be written for h{se°') and 
k{p) for k{p,ip^), etc. 

With this notation and the conformal factor 
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a conformal representation of g is defined by 

-g ^ (A^ ~g) = 2 (^^^^^-y^ y + pMp (1 " r)) / ) df (6.6) 

+(1 - fr (^^%y^ y + P- - ^)) ~ f))- 

The new coordinates do not reflect the symmetries of the underlying space- 
time, but they are sufficient to discuss the part of the space-time in the future of 
the initial hypersurface {t — 0}. We replace S by the manifold with boundary S 
introduced in section 15. 1.11 The points of dS are thought of as ideal end points 
attached to the curves p x°'{p) = pe°'{ijj^) in 5 as p ^ for fixed value of 
tjj^. The coordinates p and extend (by definition) to analytic coordinates 
on S with yO = on dS. We set 

7W' = {0 < f < 1, 0<p}, M' =M'uj+'uJ'UJ+\ 

where it is understood that the unspecified coordinate systems tp^ 'cover' the 
sphere 5^, and 

J+' ^{f = l,p> 0}, X°' = 95 = {f = 0, p = 0}, 

X' = {0 < f < l,p = 0}, X+' ={f = l,p = 0}, X'=X'UX+'. 

While the notation alludes to related sets introduced in section l531 the prime 
should warn the reader that the sets defined above differ in various aspects from 
those considered in 15.^^1 The range of p should be also bounded from above in 
these definitions. We leave this bound unspecified because its specific value is 
unimportant here, we will be concerned only with the behaviour of the metric 
in a neighbourhood of X' in M' . 

Important for the following are the observations: 

(i) the function h{se°'{tp^)) as given by the right hand side of H6.5|l and consid- 
ered as function of s and extends as a real analytic function into a domain 
where s < 0. This follows immediately from the values taken by U and its 
analyticity. 

(ii) similarly, the 1-form I given by (|6.4() extends as a real analytic function into 
a domain where p < and f > 1. This follows from 



1 \ , , 1 2{U{s) - sm)iU{s) + s 



and l|07|) with = T. 

For the following it is convenient to slightly modify the frame (|5.57|) and set 

vo = df, wi = pdp, VA = d^A, (6.7) 
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a° = df, = -dp, =dip^, A,B = 2,3. 
P 

One then gets g = cjik a* a'' with metric coefBcients 



500 = 0, 501 



.911 



Mp(i-t)) 
Hp) - 

hjpa-f)) 
Hp) 



giA 



(i-f) 



-(i-f)ph{p{i-f)) i-(i-f) 



9oA = pHpi^ - t))Ia, 

H-p{i-f)) 



Hp) 
Hp^-r)) 
Hp) 



I A, 



ffAB = {p(1-t) hip{l ~ f))y IaIb + kAsipil - f)). 

In terms of the new coordinates the metric given by H6.6|l extends analytically 
through the set J'~^ . The latter is a null hypersurface for the extended metric 
and represents future null infinity for the space-time defined by g. By contrast, 
the right hand side of Hf).6|l does not extend smoothly to 1' . However, the frame 
coefficients gik and their contravariant versions g^'^ do extend analytically to 
all of M' . It will be shown later how T' relates to (part of) the cylinder at 
space- like infinity denoted in 15. 31 bv 1. 

One has g,k = g*k + ^(p^) with 



9ik 



l+2mpf 

l + 2mpf -(l-f)(l-f f-|-4mpf2) 

fc22(0) fc23(0) 

fc32(0) fc33(0) 

so that det{g*^) < for p > 0, < f < 1 and 3^^= = g*''' + 0{p^) with 



(6.8) 



(1-t) (l-l-T-l-4 mj5f') 



1+2 mpr 





1+2 m p T 












P2(0) 
U32 



P^(0) P^(0) 






P3(0) 
33/ 



(6.9) 



Since the conformal factor A does not depend on t, the static Killing vector 
field represents a Killing field also for the metric g. In the new coordinates it 
takes the form 



i..MM{(i 



~-)df + pdp} = ph{p){{l ~ t)vo + vi}, 



(6.10) 



and extends smoothly to all oi Ai'. 

Denote by V the Levi-Civita connection of g. Since the commutators of the 
frame fields Vk vanish, the connection coefficients defined by ^tvj = VviVj = 
7i j Vk are given by the formula 



1 



7» J = 2 ^ i9ri)+ {gij ) - VI {gij ) ) 
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Again, the connection coefficients 7^ *^ j in the frame Vk extend analytically 
through {p — 0} and {f — 1}. One finds 

l^ ',^\ 9*'' {v,{gl) + v,{g*,^) - vi{g*^)) + 0{p^), 

which impUes 

7^^=T5^ {25"(,(5i,)-(l-f2)5^5i,}-f5S(5i,5S- on {p = 0}. 

(6.11) 

As a consequence of the behaviour of gij and 7^ ^ j the components of all 
tensor fields in the frame Vk which are derived by standard formulas from the 
metric and the connection coefficients, such as those of the Ricci tensor and the 
conformal Weyl tensor of g, extend analytically through J'^ and X' i.e. the 
metric g and its connection V imply in the frame Vi a smooth frame formalism 
on M'. 

It follows that the coordinate expressions of these tensor fields, such as 
Rii'ij'lg] = Rik fj.' a'' v' , and, by the argument given in 1^0] (cf. also |SHl)j 
the rescaled conformal Weyl tensor W^^ ,^'\'p'[g] = u'\'p'[g] extend 

smoothly to . Unfortunately, this does not give us the needed details about 
the components Rjk and it does not tell us anything about the behaviour of the 
frame components jki [g] of the rescaled conformal Weyl tensor on X' and the 
critical set X+ . This requires detailed calculations. Only the analyticity of h 
near i is required to control the smoothness of the fields near . This follows 
from the ellipticity of the conformal static field equations near i. To deduce the 
desired behaviour near X' , however, one will have to invoke at least, as discussed 
in section [5l)l the regularity condition (|5.89(l with p^, ^ 00. The detailed form 
of the conformal static field equations will thus become much more important. 



6.1 The Ricci tensor of g near J' 

The tensor 

is needed to integrate the conformal geodesic equations which define the setting 
introduced in section |5| The purpose of this section is to demonstrate 

Lemma 6.1 The frame components Ljk ~ L[g]pi^i j v" k extend as real an- 
alytic functions to X' with 

1 1 - f2 1 
Lok^-^S^k, Lii—t ^ — , Li^ 0, LAB^-^kAB{Q) as p — > 0. 

Proof: Under the rescaling g ^ g — A? g the tensor 
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transforms into 

L[9]p^ - L[g]p, - I Vp V. A + ^ A A gp,. 

Suppose g = dt^ + h is a. static vacuum solution and 5 = A^ g = A^ (u^ dt^ + 
h) = dt^ + h* with 

N ^Av, h* = K^h^ A^Q-^h^ /l*f,(x'=)dx°da;^ 

/i = ^(a;°), v^vix"), n^nix''), A = A{x''). 

In the foUowing the gauge H4.31|l and coordinates satisfying H4.17|l will be as- 
sumed. The connection coefficients of the metric g in the coordinates i, a;" are 
given by 

Ta c[g] = To ^ c[h*] (the Levi-Civita connection of h*), 

n ' M = 0, r/ tig] = Tt * = ^ A ^- 

and L[g\pi, is given by 

i[5]« = -^^aA^i?'^A+^Z?aAi5"A, (6.12) 

L[9\ta = L[g\at = 0, 

L[gU = -jD:DIA+^DcAD'=A hab, (6.13) 

where D and D* denote the h- and /i*-Levi-Civita connections respectively. 
With A = SIT- and the map $ defined by (j6.1(l one can determine from 
these formulas the frame coefficients 

Lik =<<P*iL[A'^g]); Vi, Vk > = 

< {L[g\tto<^)dtdt + {L[g\abO<^)dx°- dx^; v^, Vk > ■ 
With equations H4.36|l and the relations 

Dan = nia-^ Daa=il + y^)-^Daa, (6.14) 

DaDbn= ^ DaDbCJ-\.[^ ^ D^aDba, (6.15) 

(1 + y/JTay 2 V cr (1 + vfJ^r 

which are implied by 14.35|l . one gets 

TD''NDaA = 
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2s(l-2y^) (2-3v7I^)T-iD^TL)ccr 



(1 - Vm^) (1 + Va^)'^ 2(1 -77^) (1 + 771^)3 c/2 (1 + ^71^)2 j ' 



T D" ADaA = n 



2s 



1 



1 T-^D^TD^a 



T-ab 



C/2 (1 + ^71^)2 2(1 + ^/71^)3;' 



which allow us to obtain the following expressions for the Lj^. 

In the case of Lqo there occurs a cancellation of the second terms in (I6.12() , 
(|6.13|l respectively, so that (with the understanding that e° o $ = e°(?/j"^)) 



-^00 = 



(6.16) 



with 



T Da N A + vnT^/^ Dl DlAe" = 
(l-4^7rF+^cr)(s + 2) cr(l-/icr)i?afce°e'' 

(1 - V7^) (1 + V7^)4 



(1 + V7^)3 



6 cr 

" (1 + 77^)4 
(1-2 771^) 



1 D'^TDaU 



(Z7 + 



)2-l 



(1 - 01a) (1 + v7^)3 



2[/2 

2-2t/2 D^TDaU 



C/2 



C/3 



Since the term in curly brackets is of the order 0(T), the fimction Lqo extends 
smoothly to {p = 0} with Lqo as p ^ 0. 
It holds 



^01 



P 



$-(l-r)Lo 



with 



lMli = n-^r (-DaAD^A- DaN A) = 
vil 2 



1 



(1 — 3 a/TTct ) 



2/72(1 + 771^)2 (1 - 771^) (1 + 01^)'^ 



1 -2y/jra 



(1 - V7^) (1 + V7^)3 



2 D^TDaU 
IP ^ IP 
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so that Lqi extends smoothly to {p ~ 0} with ioi ^ 5 as p ^ 0. 

^0A = P (^^^'i'lA-pHl-f) {L[gUo<i>)e'^e%,, 



with 



(1 - Rab e 



1 TD'^TDaU 

T7T + 2 



C/4 t/e 



so that Lqa extends smoothly to {p — 0} with Lqa — s- as p — s- 0. 
extends smoothly to {p — 0} with in ^ — as p ^ 0. 

LlA — 

extends smoothly to {p — 0} with Lia — > as p — > 0. 

LAij = L[g]tt IaIb- {TL[g]ab) o $ e^^A e^'^^. 
extends smoothly to {p = 0} with Lyis ~^ ~j ^ab(O) as p ^ 0. 

6.2 The rescaled conformal Weyl tensor of g near X' 

In this section we shall make a few general observations concerning the rescaled 
conformal Weyl tensor and then specialize to the conformal static case. After a 
remark about the radiation field on we will analyse the smoothness of the 
rescaled conformal Weyl tensor near the set I' . 

Let ^ be a Lorentz metric and S a space-like hypersurface with unit normal 
h and induced metric h^j^i, = - n^j^fiy. We set pp,^ — h^u — fi^fiu, lv\p = 

e^!/Ap, and denote by Cup = Cpv\p[g\ fi^fi^ and c*^ = C*^^xp[g\n^' (the star 
on the right hand side indicating the dual) the n-electric and the ri-magnetic 
part of the conformal Weyl tensor respectively. The latter are symmetric, trace- 
free, and spatial, i.e. tV^ c^p = 0, h'^ c^p = 0. The conformal Weyl tensor of g is 
then given in terms of its electric and the magnetic part by (cf. |42j ') 

Cp.i,\p[g] = 2 (p^[x Cp]p. - Pp.[x ip]u - "[A c*]5 p,^ - h[p 5*]^ Ap) ■ (6.17) 

Suppose that g is a solution to the vacuum field equations. Then the first 
and second fundamental form hab and Xab induced by g on 5 satisfy the Gauss 
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and the Codazzi equation (expressing the pull-back of spatial tensors to S in 
terms of spatial coordinates x"") 

rab[h] = -Cab + Xc'' Xab " Xca Xb " , (6.18) 
DbXd(a~^c)"' = -~C*ac- (6.19) 

This allows us to express the conformal Weyl tensor in terms of hab and Xah- If 
5 is a hypersurface of time reflection symmetry, so that Xafc = 0, these equations 
imply 

rab[h] = -£ab, = 0, (6.20) 

and the Weyl tensor assumes the form 

Cy,v\p[g\ = 2 {p^^x Cp]f, ~ Pf,[x Cp]^) = -{p c)^^Ap, (6.21) 

where denotes the bi-linear Kulkarni-Nomizu product of two symmetric 2- 
tensors (cf. [5]). 

If A is an arbitrary conformal factor, the rescaled conformal Weyl tensor of 
g = A'^ g is given by v\p[g] = i/Ap[ff]. In view of the behaviour of the 

conformal Weyl tensor under conformal rescalings, one gets (observe the index 
positions) 

Wp,xp[g] ^ ACp,xp[g]. (6.22) 
Its electric part with respect to the ^-unit vector h is then given by 

WpM^^^'cpM (6.23) 

With h = il^ h, the gauge (|4.31|) . the general transformation law 

r^^lh] = rab[h] + DaDbfl + hab {^^^ DcD^n - 2 D^VL D^fl), 

and the equation 2 Ah ft = 3 Da^ D'^ft, one gets from (|6.2()|l and ()6.23(l in the 
general time reflection symmetric case 

Wab[g] = -{An)-\DaDbn-^habD,D''n + nrab[h]) on 5 = 5U{^}. 

(6.24) 

A conformal scaling which represents space-like infinity (with respect to the 
initial hypersurface S and with respect to the solution space-time) by a point 
is achieved by choosing A = f7 on S. With this particular choice one has 

Wab[g] = -i^'^Da Dbn-^habDcD^n + n rab[h]) (6.25) 

= 0(T"^/2) as T ^ unless m 0. 

We note that in the massless case the precise behaviour depends on the freely 
prescribed metric h on S near i. In the massless case one has D, — a and the 
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comparison of the expression for Wab [g] with (|4.3(j|) shows that in the case where 
h represents conformaUy static vacuum data one has 

Wabig] = -l-^rab[h], (6.26) 

i.e. the rescaled conformal Weyl tensor is smooth. 

We return to the case where A = fl T~^/^. With (|6.14ll . (|6.15() we get then 



Wabig] = T \ + VJ^) (Da DbCr - - AhCT hab) (6.27) 



i W - {3DaCr Dbcr - Dccr D" a hab) + cr (1 + JJTaf Tab 
Z V a 



- ^ {3 Da<J Dba - Dc<J a hab) - U {1 + rab[h] 
4 cr^ 2 



5- (1 + x/Jj^) ^ab, 

where we use E^h defined by the right hand side of l|4.36|) without assuming 
h to be conformaUy static. If h is conformaUy static the electric part of the 
rescaled conformal Weyl tensor on S is given by the right hand side of H6.27|l 
with Eafe = 0. In the present conformal gauge, defined by (|4.31() . one has 

Sa6 = 0(T3/2) near i, 

for any time reflection symmetric initial data h. 

If the solution is static and written again in the form g = v'^ dt^ + h, then 
equations (|6.17(l , H6.20|l hold with h ^ ^ dt and ^-independent fields for each 
slice S = {t = t^} with = const. The relations above then imply for all 

{t,x°') 

W^,xp[g] = -T-i(p w)f,,xp (6.28) 

with 

Pfiiy — hf^i/ Tip Til/, — ^ 

With (g) denoting the tensor product, we write for arbitrary 1-forms a, c 
a(^sC — a(^c + ci^a, = a ® a, 
and note that the Kulkarni-Nomizu product is symmetric, i.e. 

m(2)n = n(2)m, (6.29) 
for symmetric 2-tensors m, n, and satisfies for arbitrary 1-forms a, c, e 

(a (g) a) (a (g)s c) = 0, (a 0s e) (a 0s c) = -(a a) (c 0s e). (6.30) 

We show how it follows in the present setting that the radiation field van- 
ishes on J'^ . Since the extended Killing vector field K is tangent to the null 
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generators of without vanishing there, the complete information on the 
radiation field is contained in the field 



= -J-^K" KPp^pWabdx^dx^ = -T-^ {vftf Wabdx" dx'' 
^ m (1 - Vga)^ \^{2shab-3a-' DaaDha)dx''dx'' 




+ ^(1 + rat[h] dx" dx^ 



with s as given in l|4.32|l . Because of the relation 



a-^DaaDb(7 = A U'^ {U^ e, et - U T^/^ (e, D^U + eb DaU) + T DaU Db U), 



and the factor a in the second term it follows that 

K"" KPWp^Xp[g]dxPdx^ ^ -2mp^df^ as f 1, p > 0. (6.32) 

Thus, the pull-back of K"^ Rp Wfj.,y\p[g] dx^ dx^ to which provides the ra- 
diation field up to a scaling, vanishes everywhere on 

Lemma 6.2 The components Wijki[g\ — f^^^ Cpiy'y p'[g]vP- iv" j kvP i of 
the resettled eonformal Weyl tensor of g in the frame Vk extend as analytic 
functions to X' . 

Proof: In the coordinates x^ given by H6.1|l the rescaled eonformal Weyl 
tensor is obtained as the product of 



(6.31) 



(To$)-i 



(P(l-r)) 



-2 



with the Nomizu-Kulkarni product of 



pj = (hab 0^)dx°' dx^ -{(y^)o ^fdt^ ^ p'^ + P2 + P3 + P4> 



and w' — w[ + W2 + w'^ + w'^ + w'^, where 



p[ = -2{il -f) dp- pdf)\ 





Pi = P' 



f2(l-f)2fc, 
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and 



' 1 ' I ' I ' 

Wi + Wn + Wo + Wy, 




with 



^({^(.+6f,-,}.*)ri. 
:^-({f^/=.}o*).. 



( 



{ 



3 m 
Wa 



I o$ ) ((1 -f)dp-pdf) 



We used above the relation (|6.31() and set 



The desired result on the behaviour of the rescaled conformal Weyl tensor 
near I' is obtained now by showing that for arbitrary frame vector fields w„ one 
has 

<p' (Z) w'; v„ Vj, Vk, VI, > = 0{p^ (1 - f)2). 
From it follows that 

p'l (Z)w[= p[ = p'2 0w[ = 0. 

Observing H4.27(l one finds by inspection 

<p[]v,,Vj>^0{p^), <PM;v^,Vj>^0{p^{l-ff), for M = 2, 3, 4, 

< -u;4; -y,, >= 0(p2 (1 - f)2), < w'j^; v^, Vj > ^ 0{1), forA^ = 2,3,5, 



<p'M^w'j^;v^,Vj,Vk,vi>^0{p^{l-f f) for M = 2, 3, 4, TV = 2, 3, 4, 5, 

< p[ w'^; v„ Vj, Vk, VI >^ 0{p^ (1 - ff). 
The remaining term is given by 



and thus 



p'l (w^ + w^;) + (p^ + p'^) w[ ^p'^Cdm 



with 




^ ({[/ (3 s [/2 + 6)} o $) fc - ({^ [/ (1 + VTI^)' o $) dx'^dx^ 
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4 > ' \p{i~T){vn){p) ^ p(i-t) ; 

For the three summands to be considered here we get the following. From 
3 s [/2 + 6 = 0(T) it follows that 

< P'l {{U (3 s + 6)} o $) fc; wj, t-fc, t-z > = 0(p2 (i _ f)^). 

Because of 

-ip'i e-^ + p (1 - f) e'"^ de'') ®, ((1 -T)dp- pdr), 
it follows by that 

<p'i ({y (1 + VJj^)'^ rab[h]^ o <I>^ dx^dx^] V,, Wj, Wfc, -w; > = 

0(p2 (l-ff). 

It holds that sU^ + QU ~ 0{\) and by inspection it follows that 

\p{i-T){vn){p) / 

= 0(pMl-f)2). 

7 Static vacuum solutions near the cylinder at 
space-like infinity 

The conformal extension considered in the previous section relies on specific 
features of static fields. We use it to show that the construction of the cylinder 
at space-like infinity in section[Sl which is based on general concepts and applies 
to general solutions, is for static vacuum solutions as smooth as can be expected. 



Theorem 7.1 For static vacuum solutions which are asymptotically flat the 
construction of section \^ is analytic in the sense that in the frame j5.57| j all 
conformal fields, including the rescaled conformal Weyl tensor, extend to ana- 
lytic fields on some neighbourhood O ofX in Af. This statement does not depend 
on a particular choice of (analytic) scaling of the (analytic) free datum h on S. 

This result will be obtained as a consequence of Lemmas 17.21 [731 and 17. 41 below. 

The construction of section [S] will be discussed here for static solutions in 
terms of the initial data h and in the gauge given by (|4.31|) . and the field g 
given on A^' in the coordinates defined by H6.1|l . The effect of a rescaling of h 
will be dicussed seperately because it is of interest in itself. 
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The conformal factor 8 is assumed in the form (|5.44ll . 15.46|l with 



K = uj = 2Q \DanD''n\-i = 2n{i + ^/jr^y^ \/2\s\a. (7.1) 

It foUows that uj = + 0(T) and 9 = 1/2 \Da^D''n\^ = T^/^ + 0(T) so 
that 

on 5: hm T^/^cj^i = hmcjA-i = hm eA-i 1. (7.2) 

The metric h induced by g = 0^ g on 5 is given hy lo~^ h. 

The main ingredient of the gauge for the evolution equations used in section 
I5.2l are the conformal geodesies generating the conformal Gauss system described 
in section We shall try to control their evolution on tW' near X' . Following 
the prescription in section [2. II we assume that the tangent vectors x — dx/dr 
of the conformal geodesies with parameter r satisfy 

x±S, e^g{x,x) = l on S. (7.3) 

With the frame (|6.7|) and the coordinates (|6.1|l this translates into the initial 
condition 

X = ^ df + u{p) dp = ^ (vo + vi) = X'v, at f = 0, (7.4) 
P P 

with 

X' = 5\ + 5\ + 0{p) as p -> 0. 
For the following we need to observe besides g ~ h? g the relations 

3 = ^3 = 9^5, with HeeA^^O. 

For the connections V, V, and V of g, and g respectively we have relations 

V = V + 5(/), V = V + 5(/), V = V + 5(/), 

V = V + S'(e-ide),, V = V + 5'(A-idA). 

The comparison gives f = f — dQ and f =^ f — A^^ dA which imply the 
relation 

f^f + e-^de-A-^dA = f + n-Un, (7.5) 

between the 1-form / which is obtained if the conformal geodesic equations are 
written in terms of the metric g, the 1-form / which is supplied by the conformal 
geodesic equations written in terms of the metric g, and the conformal factor 
which relates g to g. 

By the choices of section we have < /, i > = everywhere on the space- 
time and <d8,i>=0on5. Since A has been chosen to be independent of t 
and dt is orthogonal S, it follows that < dA,x >= and thus < /, i > = on 
iS*. Observing the pull back of / to 5 given by (|5.45|l and the relation 

n = T^^^uj~^ on 5, (7.6) 



76 



we find that the pull back of / to 5 is given by 1/2 T ^ dT. From this one gets 
in the frame (|6.7I) and the coordinates (|6.1II with f = 

f = {1/2 T-^DaT)o^ dx'' ^ f^a' with /, = -5" , + on 5. (7.7) 

The relation < /, i > = and equation (|7.5|l imply the ODE 

ri = n</,i>, (7.8) 

along the conformal geodesies, which, together with H7.6|l . will allow one to 
determine 11 once < /,i > is known. 

7.1 The extended conformal geodesic equation on / 

With respect to the metric (|6.6|l a solution to the conformal geodesic equations 
is given by a space-time curve x'^{t) = (f (r), /o(t), %P^{t)) and along that curve 
a vector field X{t) and a 1-form /(r) such that 

x = X, 

VxX = -2 <f,X>X + g{X, X) /», 

Vx / = < /, X > / - i g{f, J) X' + L{X, . ). 

With the expansions X = Vi, / = fia^, g = gjkCt^a'', = g-''^VjVk, 
L = Ljk , = k Vi, the equations above take in the domain where p > 
the form 

^f = xo, ^-p=-px\ ^^^^x\ 

dr dr dr 

which is the equation x^ = X^ i, relating the coordinate to the frame expres- 
sions, 

+ 7/ fe X^ X'' = -2fkX^X'+ g,k X^ X^ f fi , 

a T 

-^fk-lj\ e h^IiX'h~\ fi h 9ik X'+L,k X^ . 

CLT Z 

Note that the functions cjjk^ g*', 7j*i, Ljk entering these equations extend 
by analyticity through T' into a domain where p < 0. Assuming such an exten- 
sion, we get the extended conformal geodesic equations. Since also the data are 
analytic on S, it makes sense to consider these equations in a neighbourhood of 
J'. 

Lemma 7.2 With the values of Ljk on J.' found in Lemma \6.1V the initial data 
X = (0,0, V^') and (cf {TJ); X* = 5\ + 5' f, = ~5\ + 5\ on 

1° determine a solution x{t), X{t), f{T) of the extended conformal geodesic 
equations with t = on 2'^ and 
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By analyticity it extends as a solution into a domain < t < 1 + 2e for some 
e > 0. The extension to T' of the conformal factor 11 which is determined by 
1 7. 6] ) and i7.8j) takes the value 11 = 1 on T' . 

Proof: With the ansatz x(r) = (f(r), 0, V--^'), X{t) ^ X°{t) vq + X^{t) vi, 
f = /o(r)a*' + fi{T)a^ those of the extended conformal geodesic equations 
which are not identically satisfied because of H6.8|l . H6.9|l . H6.11|l are given by 

dr 

X" + 2f X° X^ - f {1 ~ f^) X^ X^ 

dr 

= -2 (/o X° + h X') X" + (2 X° X'-{1- f')X' X') ((1 - f 2) /o + A), 



— X^~fX^X^ = -2 (fo X° + fi X^) X^ + (2X°X^ - (I - f^)X^ X^) fo, 
dr 



-^fo-ffoX'^ ifo X" + h X') fo-l ((1 - f') fofo + 2 /o A) X' + ^ X\ 



^fi-f fo (XO - (1 - f 2) X') + f h X^ = (/o + h X') h 

^\{{l~f^)kh + 2hh){X^-{l-f^)X^) + \x'-^—:^X\ 

A calculation shows that the solution of this system for the prescribed initial 
is given by 

f = r, X0 = 1, X' = -^, /o = --i-^, /-=1. (7.9) 

This proves the first assertion. With the solution above equation (|7.8|) reads 
n = and we have 11 = 1 on /° by H7.2|l . This proves the second assertion. 

Remark: The ODE above is sufficiently complicated so that giving the solution 
explicitly deserves an explanation. In 15.62|l is given the conformal factor and 
the conformal representation of Minkowski space which result from the general 
procedure of section In H5.63|l is given the coordinate transformation which, 
together with the conformal factor, relates the conformal metric to the standard 
representation of Minkowski space in coordinates t and r. 

If the Minkowski values to = 0, C/ = 1, hab = —Sab are assumed in section 
the metric g reduces by H6.8|l to the metric g*^ a* with m = 0. One can 
consider this as the lowest order (in p) approximation of the general version of 
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g. Tracing back how the functions r, p, A in section El are related in the fiat 
case to t and r, one finds 

1 ^ . 1 

r — —, rr, t — —, rr, A = -. 

p(l-r) p(1-t) r 

The conformal factors in the conformal representations thus agree but the co- 
ordinates are related by the transformation 



P = p(1 + t). (7.10) 



This implies 



P \ P J P \ P J 

(7.11) 

The left hand side is the conformal Minkowski metric (|6.8|l with m = while 
the right hand side is the metric g* given by H5.62|l . The conformal geodesies 
underlying H5.62|l have tangent vector X = dr and 1-form / = With H7.10|l 
these transform into 

X ^ df + . ^ _ pdp = vo + -^-^ vi, 
1 + T 1 + r 

/ = — - -r— ^flT" = --r— ao + ^i, 
p I + T 1 + r 

from which one can read off 1)7.9(1 . 

The hypersurfaces {p — p^ = const. > 0} are in general time-like for the 
metric The form of g^ suggests that these hypersurfaces approximate null 

hypersurfaces in the limit as p^^^ ^ 0, but the conclusion is delicate because of 
the degeneracy of g'^ on X' . The discussion above shows that they do become 
null asymptotically in the sense that for the metric on the left hand side of H7.11|l 
the hypersurfaces {p = const.} are in fact null. To some extent this explains 
why the coordinates given by 1(6. l|l had a chance to extend smoothly to J'^ and 
to provide a description of the cylinder at space-like infinity. 



7.2 The smoothness of the gauge of section [31 for static 
asymptotically flat vacuum solution near I 

Let Sext denote an analytic extension of S into a range where p < so that 
p, ip^ extend to analytic coordinates. If the set Sext \ 5 is sufficiently small, 
the following statements make sense. The initial conditions ((7.4(1 . 1(7.7(1 extend 
analytically to Sext and determine near Sext an analytic congruence of solutions 
to the extended conformal geodesic equations. It therefore follows from Lemma 
17.21 and well known results on ODE's that, with the e of Lemma 17.21 there 
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exists a p# > such that for initial data r(0) = 0, p(0) = p' with \p'\ < p^, 
V'^(O) = ifi"^ and those implied at these points by H7.4|l . (|7.7|) the solution 

f = f(T,p',i;^'), p = p{T,p\ij^'), V^^ = V^^(r,p',V^'), 

- X\t,p',^^'), h = /fe(T,p',V'^'), 

of the extended conformal geodesic equations exists for the values < r < 1 + e 
of their natural parameter and the function 11 is positive in the given range of 
p' and T. 

Taking a derivative of the equation satisfied by p and observing H7.9|l gives 
d (dp \_fdp. \ 1 

which implies by H7.9fl 

It follows that the Jacobian of the analytic map 

{t,p',^^')^x'^{t,p',^^'), 

takes the value 1 + r on I' and for sufficiently small p# > the Jacobian does 
not vanish in the range < r < 1 + e, \p\ < p^. The relations A = 011^^, 
n > 0, and 8 — {uj~^ il)* (l — r^) imply that the curves with p' > cross J'^ 
for T = 1. It follows that t, p' , and define an analytic coordinate system 
in a certain neighbourhood O' of I' in Ai' , such that (suppressing again the 
upper bounds for p') O' J+' = {r 1, p' > 0}, X' = {0 < r < 1, p' = 0}, 
X+ = {t = 1, p' = 0}, and C is ruled by conformal geodesies. 

The metric g — g, the connection coefficients of the connection V and 
the tensor fields (cf. H2.11|l ) 

Lf_iu — [g] — fu + ffi fi^ — 2 dtif f\ 1^1 

/ = / - n- 1 d n, ly^.pA [5] = n w^^.p^ [g] . 

in the frame (|6.7() extend in the new coordinates as analytic fields to O' . 

Given these structures and the conformal geodesies on C, the construction 
of the manifold N as decribed in section [3 poses no problems. With the given 
analytic initial data on S it only involves solving linear ODE's corresponding to 
(|2.21|l . such as 

e% + 7/ / X' e' k = -Ii X'e\- fi e' X' + g.i X^ e' k 5^™ /m, (7.12) 
or its spinor analogue, along the conformal geodesies. This allows us to conclude 

Lemma 7.3 Starting with static asymptotically flat initial data in the gauge 
\4-31\l , the construction of section\^ leads to a conformal representation of the 
static vacuum space-time which is real analytic in a neighbourhood O of the set 
I inM. 
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7.3 Changing the conformal gauge on the initial slice 

It will be shown now how the construction described in section [S] depends for 
static vacuum solutions on rescalings 



UJ 



h ^ h' = uj-^ h, n^n'=^dn on s, 



with analytic, positive conformal factors i?. 

There are harmless consequences such as the change of the normal coordi- 
nates x'^ ^ X = X "■{x'^) with X ^{0) — and a related change Cq ^ = 
•d'^ s"^ a Cc of the frame vector fields tangent to S. Here s"^ a denotes an analytic 
function on S with values in 50(3) such that s"^ q ^ J'^ a as p ^ 0. These 
changes will simply be propagated along the new conformal geodesies. 

Critical is the transition from the congruence of conformal geodesies related 
to f2 (the Q- congruence) to the new one related to f2' (the ^' -congruence) . If the 
curves are considered as point sets, the two families of curves will be different 

if n'-^ dn' ^ n-^ dn ^ d-^ d-d ^ o {d. gO]). 

The rescaling above implies on S the transitions 

\\dn\\h ^ Wdn'Wh' =^\\dn\\h, 

2n , 2n' u;S 

UJ = 



Is ^ 



with the function 



which extends to 5 as a analytic function of p and 

It follows from the initial conditions for the f2-congruence that 

C^x±S, Q'^g{C^x,r^x) = 1, (7.13) 

and for the transformed 1-form that 

< /', i > = 0, fs^f's^ uj'-^ dJ = fs + ^-^d§- r ' di, 

where the subscripts indicate the pull back to S. These two lines give the initial 
data for the fi'-congruence if the conformal geodesic equations are expressed 
with respect to the rescaled metric g' and its connection V'. 

To compare the ri'-congruence with the 0-congruence we observe the con- 
formal invariance of conformal geodesies (cf. ^j) and express the equations for 
the i7'-congruence in terms oi g and its connection V. The space-time curves, in- 
cluding their parameter r', then remain unchanged. The 1-form is transformed 
because of g = (6 Q'~^)^ g' according to f'^f*^f- (6 Q''^)-^ d (9 G'-^), 
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which imphes < f* ,x > — 0, = fg + d ^ d-d on S. If this 1-form is expressed 
in terms of the g-orthonormal frame with cq J- 5, one finds 

/*=</*, eo> = 0, /:=< r,ea> = fa + ^-' <di9,ea>, a=l,2,3. 

(7.14) 

The fields ^ ^ x, are the initial data for the rj'-congruence in terms oi g, e^, 
and V. Since C ^ 1 and < d-d,ea > = 0{p) as p ^ 0, it follows that 

l^^l, i-r'i^O, /*-/^->0 as p^O. 

As a consequence, the initial data for the Q!- and the fi-congruence have coin- 
ciding limits on and the corresponding curves are identical on X'. 

Assuming now the conditions of section [Y. 21 and using arguments similar to 
the ones used there, we conclude that in a certain neighbourhood O' of X' in 
M.' the gauge related to the ll'-congruence is as smooth and regular as the one 
related to the fi-congruence. Thus we have 

Lemma 7.4 In the case of statie asymptotically flat space-times the construc- 
tion of the set X' is independent of the choice of Q and the set X' introduced in 
section\^ coincides with the projection 7r'(X) of the cylinder at space-like infinity 
as defined in section\^ 

We note that the comparison of the 17'- with the il-congruence leads in the 
case where the solution is not static and thus not necessarily analytic still to 
similar results if the solution acquires a certain smoothness near J'^ U X* . In 
the case of low smoothness, however, the detailed behaviour of the different 
congruences needs to be analysed in the context of an existence theorem. 



8 Concluding remarks 

Concerning the regularity conditions we have now the following situation. For 
static asymptotically flat solutions with m ^ the conformal extensions to X 
are smooth (in the sense discussed above) and their data satisfy the regularity 
condition H5.89|l withp* = oo. In the massless case condition H5.89|l withp* = oo 
is necessary and sufficient for space-like infinity to be represented by a regular 
point in a smooth conformal extension. In the general time reflection case with 
m 7^ conditions (|5.89() are necessary but not sufficient for the s-jets Jjiu), 
p e N, to be regular at the critical sets X*. Thus, the mass m — 2W{i) and 
also the derivatives of d"aW{i), a gN'^, play a crucial role for the behaviour of 
the Jjiu) at X^. The mechanism which decides on the smoothness remains to 
be understood. 

Only the d-jet Jjo{u) and the s-jets J^~^{u) are needed to obtain Jj{u) 
by integrating the transport equations on X. Since the left hand sides of the 
transport equations are universal in the sense that they do not depend on the 
data, it follows that (it) is uniquely determined by Jjq{u) for p G N. In the 
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static case the s-jets Jjiu) are regular. In jJS] has been exhibited a class of 
data which are asymptotically static of order p for given p S NU {00} and which 
are essentially arbitrary on given compact sets. It follows that for prescribed 
differentiability order p there exists a large class of data for which the s-jet Jjiu) 
is regular on X. 

We expect there to be a threshold in p beyond which the regularity of Jj{u) 
ensures peeling resp. asymptotic smoothness of a given order of differentiability 
and below which the singularity of Jj{u) implies a failure of peeling. This order 
is likely to be low enough such that the behaviour of J|(u) with q < p can 
be controlled by a direct, though tedious, calculation. However, if asymptotic 
staticity does play a role here, one should try to understand the underlying 
mechanism. It would be quite a remarkable feature of Einstein's equations if 
asymptotic staticity could be deduced from asymptotic regularity at null infinity. 
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